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§ 1. Object and Aims of the Paper. 

The usual solution for the extension and compression of elastic bars assumes that the 
latter are strained under a normal tension or pressure uniformly distributed across 
the plane ends. In like manner the solution for torsion of such bars assumes that 
the external forces which cause the torsion consist of a determinate system of 
tangential stresses, acting across the plane ends. 

In both cases the solution is such that the torsion and extension are transmitted 
throughout the bar without change of type. Such terminal conditions of stress, 
however, do not usually occur in practice, and it accordingly becomes of considerable 
interest to find out how the results obtained for such a theoretical system of loading 
are modified, if at all, when we consider applied external stresses which give a closer 
representation of every-day mechanical conditions. 

The present paper is an attempt towards the solution of this problem in three 
cases, which appear of especial practical interest. 

The first case is that of a bar which is subjected to a determinate system, of normal 
radial pressures and of axial shears all over the curved surface, the radial pressures 
being symmetrical about the mid-section and the shears having their sign changed. 
Thus the cylinder is subjected to a total axial pull, due to the shears, and also to a 
given transverse pressure. The plane ends are free from stress, except for a self- 
equilibrating system of radial shears, which will have little or no effect at points at 
some distance from the ends. 

A special case worked out is that where the normal pressure is zero throughout, 
but a determinate axial shear, which has been taken constant, is made to act over 
two equal rings on the surface of the cylinder. 

This will give us valuable information about a system of stress which often occurs 
in practice, in testing machines, for example, in which a specimen is pulled apart by 
means of pressures applied to the inner rims of projecting collars (see fig. 1). The 
shaded parts of the figure represent the "grips," and if S be the total pull, applied, 
this is transmitted to the test piece by means of pressure applied along CA, C'A . 
Now consider the thinner cylinder in the middle ideally produced inside the thicker 
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ends. It is in equilibrium under the stresses, radial and tangential, between the 
inner core and the hollow cylinder produced by the revolution of ABDC. 

Fig. 1. 







But what are these radial and tangential stresses ? If we consider the equilibrium 
of the outer hollow cylinder only, we see that the resultant of the stresses across 
AB, A'B' must exactly balance the pull S, however applied. The radial stress will 
probably be small, as it has no external traction to balance, and the longitudinal 
shears are therefore equivalent to S. Thus the thin cylinder inside is really 
stretched, not by normal traction over the flat ends, but by longitudinal shears over 
the curved surface, and a careful investigation will show that, in every practical case, 
extension is obtained by the application of an axial shear to the curved surface of the 
cylinder, never of tractions to the flat ends. The general effects of such a distribution 
appear, therefore, of great practical interest. 

The second problem discussed is that of a cylinder of moderate length, which is 
compressed between two rough rigid planes in such a way that the terminal cross- 
sections are constrained to remain plane, but are not allowed to expand, their 
perimeter being kept fixed. By adding a suitable uniform distribution of pressure to 
a load system of this type, we can obtain the solution for a cylinder constrained in 
such a way that its ends expand by a definite amount. These two problems are of 
importance with reference to the behaviour of a block of stone or masonry when 
tested between millboard or metal planes, which practically hinder the block from 
expanding, and when tested between sheets of lead, which, on the other hand, favour 
the expansion of the block. The widely divergent results obtained for the strength 
of the same material when tested by these two methods have troubled many 
elasticians. Unwin ('Testing of Materials of Construction') is of opinion that the 
results obtained when sheet lead is used are unreliable ; whereas Professor Perry, in 
his ' Applied Mechanics,' states that the true strength of the materia] is the one 
given by the lead experiments, and should be usually taken as half the published 
strength. 

Finally, the third problem treated is that of a cylinder subjected to transverse 
shears over the parts of the curved surface near the ends, these shears being 
equivalent to a torsion couple. This is really the analogue, for torsion, of the first 
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problem for tension, and corresponds to the case of a bar "gripped" as explained 
above, and twisted. This again is the method by which torsion is practically 
produced in most cases- — almost always in laboratory experiments. 

5 2. Method of Solution Adopted* Historical References, 

The method adopted has been to obtain symmetrical solutions of the equations ot 
elasticity in cylindrical co-ordinates and to express the typical term in the form 

C ? 8 \kz] Xf(r) 9 

r, <-/>, z being the usual cylindrical co-ordinates. 

The expressions for the strains and stresses, over any coaxial cylinder, are therefore 
series of sines and cosines of multiples of z. The arbitrary constants of the coefficients 
are determined by comparison with the coefficients of the Fourier's series which 
express the applied stresses at the external boundary. 

This method is not a new one. It has been indicated by Lame and Clapeyrox 
(" Memoire sur fequilibre interieur des corps solides homogenes," c Crelle's Journal/ 
vol. 7), but it has been for the first time worked out with any completeness by 
Professor L. Pochhammer (" Beitrag zur Theorie cler Biegung des Kreiscylinders/' 
c Crelle's Journal/ vol. 81, 1876). Professor Pochhammer obtains the general 
solution of the elastic equations for an infinite circular cylinder subject to any system 
of surface loading, repeated at regular intervals. This he applies to the case of a 
built-in beam. The solution is not restricted to be symmetrical about the axis of the 
cylinder, but is perfectly general. The complete accurate expressions are, however, 
quite unwieldy ; but, as the result of expanding the functions involved to the first 
two or three terms, Professor Pochhammer obtains far more manageable expressions, 
which he is eventually able to identify with those previously given by Navier and 
de Saint Venaxt for more special cases of loading. It is to be noted, however, that 
Pochhammer restricts himself solely to the case of bending, and that his approxima- 
tions depend upon the ratio of diameter to length being a small quantity. 

The same general expressions have been independently arrived at by Mr. C Chree 
( Ci The Equations of an Isotropic Elastic Solid in Polar and Cylindrical Co-ordinates, 
their Solution and Application," ' Camb. Phil. Trans./ vol. 14). Here, again, the 
solutions are not restricted to be symmetrical. The symmetrical terms, however, 
agree with the solutions of the present paper, but the latter are obtained by a process 
slightly different from that of Mr. Chree. Mr. Chree has also given a solution ol 
the symmetrical case proceeding in powers of r and z. Using each form of solution 
independently, it is not possible to satisfy the condition that there shall be no stress 
at all on the curved surface ; this is effected in the second problem of this paper, by 
means of a combination of the two types of solution. 

In the paper referred to, Mr. Chree, like Professor Pochhammer, has not, so far 
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as I am aware, applied his general solution to the problems of tension and compression. 
He does give one example of torsion, which he obtains by applying an arbitrary 
system of cross-radial shears across the flat ends. Such a system, we have seen, 
would not usually correspond to what occurs in practice. 

Mr. Celree has written several other papers (" On some Compound Vibrating 
Systems/' ' Camb, Phil. Trans,/ vol. 15, Part II. ; "On Longitudinal Vibrations/ 5 
'Quarterly Journal of Mathematics/ 1889; " Longitudinal Vibrations in Solid and 
Hollow Cylinders," 'Phil. Mag./ 1899; "On Long Rotating Circular Cylinders," 
'Camb. Phil. Soc. Proc./ vol. 7, Part VI., &c), which deal with the solutions of the 
equations of elasticity in cylindrical co-ordinates, with special application to vibra- 
tions and rotating shafts ; but I cannot find that he has anywhere returned to the 
statical problem and its solution by means of sine and cosine expansions. 

[October 3, 1901.— Professor Schife ('Journal de Liouville/ Serie 3, vol. 9, 1883) 
has attempted the solution of the problem of the cylinder compressed between 
parallel planes, which is one of those treated of in the present paper. His solution 
is expressed in a series, not of circular functions, but of hyperbolic sines and cosines 
of nz, the successive values of n being obtained as roots of a certain transcendental 
equation. This enables him to satisfy the conditions at the curved surface, but the 
arbitrary coefficients are finally determined by the conditions over the plane ends. 
He assumes both the radial shear and the molecular rotation in a diametral plane 
to be given by known functions, f(r) and F(r), over the plane ends, and from these 
he succeeds in obtaining the coefficients. As he has only a single set of the latter 
left to carry out the identification, his functions f(r) and F(r) are not really inde- 
pendent. Theoretically only the shear f(r) should be required, and in a practical 
problem even this is unknown, the total pressure being all that is given. The actual 
distribution of this pressure does not appear to enter into Professor Schiff's solution. 
Also the fact that the values of n are roots of a transcendental equation singularly 
complicates the solution from a numerical point of view, and Professor Schiff appears 
to have made no attempt to translate his results into numbers.]^ 

It has therefore appeared worth while to apply the solutions involving circular 
functions of z to problems such as those sketched above. 

Of each of these I have given a concrete numerical example. Indeed, the greater 
part of the work has been spent on these numerical examples. The labour of calcu- 
lation has in most cases been considerable, owing to the slow convergence of many of 
the series involved, which has necessitated special methods of approximation. 

* Since writing the above, I find that the problem of the circular cylinder under a symmetrical strain 
has been considered by J. Thomae in two papers ( " Uber eine einfache Aufgabe aus der Theorie cler 
Elastic! tiit," ' Leipzig Eerichte,' vols. 37-38). The author has used expansions in sines and cosines of kz, 
but, as far as I can make out, the only problem he considers is that of the vertical pillar under its own 
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3. General Solution for a Symmetrical Strain, 



Let v, </>, z be the usual cylindrical co-ordinates ; also, following the notation of Tod- 
hunter and Pearson's ' History of Elasticity/ let st denote the stress, parallel to ds 9 
across an element of surface perpendicular to dt y s, i standing for any two of the 
letters r, </>, z. 

Let U) v, w denote the radial, cross-radial, and longitudinal displacements respec- 
tively, then we have (Lame, ' Legons sur TElasticite ? ), if u> v, w are independent 
of <j> : 
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d i v> 



dhi 
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X and /x being the elastic constants of Lame. 

We see from the above that <j>z and v<f> depend only on i\ the other stresses only 
on u and w. Also the equation (2) contains v only, (l) and (3) contain u and w only. 
The solution for transverse displacements is therefore absolutely independent of the 
solution for radial and longitudinal displacements. 

7 "1 7 7 

Let us now denote the operators 7 - r by £ 3 , and -- by D. Differentiate (1) 

with regard to z and (3) with regard to r 9 and remember that the order of the symbols 
I) and $■ is indifferent, then we find 



(#> -f D a ) a = . 



dn 
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(6). 
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Eliminate either du/dz or dw/dr between (6) and (7) : it is found that either of 
these quantities satisfies the partial differential equation 

(£ 3 + B^fy^O (8). 

The whole problem of the determination of the elastical equilibrium of the circular 
cylinder under any symmetrical system of stress depends therefore on the solution of 
this differential equation. 

§ 4. Solution of the Differential Equation. 

The differential equation 

(£ 3 + D 3 ) y = 

is really identical with Laplace's equation in cylindrical co-ordinates, namely, 

1 d rdY 1 (PY d*Y 

r dr dr r 2 d(f> 9j dz 2 

Suppose V independent of </> and differentiate with regard to r, we have 

(#> + D 3 ) dY/dr = 0. 

If therefore V be such a solution of Laplace's equation, y = dY/dr will be a solution 
of the given differential equation. For our purpose, however, it will be simpler to 
proceed from the equation itself. 
Assume a typical solution 

y 1 = E 1 .Z 1 , 

where R x is a function of r only, Z x a function of z only. 
We find easily 

73 i & t-^x ~— : U , ■ t . . . . . ' + . (10). 

The solutions of (9) are of the form 

I x (kr) and K x (Ar), 
where 

CO ,- t M + 2s 

t / r \ v _____ __a_.._ ._.„ 

x »v / 7 2» +2 *n («) n (s + w) J 

K/ \ / -v 1.3, (2n — 1) f^cos (x sinh <f>) ., . 

(See Gray and Mathews, 'Bessel's Functions/ pp. 66-7.) 
vol. oxcviii. — a, x 
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Consider now the equation 

(#> + J)*fy - o. 

Let y i be any solution of the equation 

(P. + D 3 ) ?/ = 0. 
Then if y be a solution of 

it is also a solution of 

Now if ^ = Ii 1 Z 1 what is the condition that we can obtain a second product solution 
y, = R,Z, ? 

We have, substituting 

A I . 3 __i,._J /11\ 

p I rr p y ' ' ' ' • ' ' ' • \ 1 /) 

Al rJ Al At 

or function of r only + function of z only = product function in r and z. 

If (11) is to be identically satisfied, this product function must be a function of r 
only or of z only. 

Case (i.). Z 3 = aZ i where a is a constant. 

We find 

or, since R x is a solution of 

which is the same as (9), R 2 is a solution of 

(3* _ p)*R a = 0, 

which is not at the same time a solution of (13). 
Now the solutions of this equation are 



But 



and similarly 



I\ (kr), K x (kr), ™ I, (£r), ^ K 2 (Ar). 



^ I x (&r) = rl/ (£r) = rl (kr) - ~ I x (/c?\) 



elk 



4 K i (^ r ) = rK o ( kr ) - j: K \ ( kr ) • 



The four independent integrals are therefore 

I x (&r) , K x (£r) , rl (£r) , ?'K (fcr) , 
and therefore the required values of ~R. 2 are rl , rK . 
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Case (ii.). R 2 = bR x ; we find, using (13), 

and therefore Z 2 is a solution of (D 3 -f /c 2 ) 2 Z 2 = 0, which is not at the same time a 

solution of 

(D* + F)Z 3 = 0. 

The possible values of Z 2 are 2; cos &2, z sin yb. 

Hence the possible sets of product functions satisfying the equation 



are as follows : 



y = A cos (kz + a) Ii (kv) 
B cos (/b + /3) Kj (/it) 
C cos (yfe -f- y) rl (kr) 
D cos (yb + S) r.K (At) 
E?j cos (fe: + e) Ii (&>') 
Fz cos (/b + 0) K x (&r) 



(14). 



§5. Solution under given conditions of Surface -loading ; the first problem. 

Let us now consider first the case of a circular cylinder under the following system 

of stress : 

rr/fi = a given even function of z (= f(z) ) over the curved surface v = a , 
r£//x = a given odd function of 2: ( = 1// (z) ) over the curved surface r = a 9 
zz = over the plane ends ^=ic 

Since du/dz, divjdr are both solutions of (8) we may have them composed of a series 
of terms as follows : 



du __ v f A x cos (fc + a x ) l^kr) + C x cos (&s + y } ) rl (hr) 
dz ~l ' + E^ cos (&s + Cl ) I x (£r) 



► . . . (15). 



(lw = ^ f A 2 cos (As + a 2 ) I^/tT) + C 2 cos (/te + y 2 ) rI (Ar) 1 /-^ 

*' I + E 2 z cos (fe + €3) Ii (£r) J ' l ; ' 

No K-functions have been introduced in this case, as they lead to infinite terms at 
the axis. 

x 2 
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T? 



Also the conditions of the problem require that u shall be an even function of % 
and to an odd function of z. Hence 



a 



etc 



2 



Tl = 7; 



7T/2 , €, 



€o 



0. 



Integrating (15) and (16) we have 






f A. , E / cos 7## N 

X (?*) + S \ — - a cos Znzl-, (lev) — y 1 cos kz,rl (kr) + ~ 2; sin &z -f- — '— 



Ii (*r) > 



w = (2) + S<J ^ 3 sin fel (&r) + y sin &z. rlj (&r) + -j % cos /b I (kr) > 



IV 



To find the relations between the constants we must substitute in equations (l) and 
(3). We then find the following relations : 



&\ (r) = , 



• (17). 



DW (z) = . 



(18) : 



(A x - A g ) (X + /*) /c 3 + 2/>;{C 1 (X + 2fi) - /uE) - (X + fi) EJ = . . (19), 
(Aj - A 2 ) (X + /*) £ 3 + 2Jk{C 1 (\ + fi) + /xCg - (X + 2/x) E 3 } = . . (20), 



C x = C 2 = C, say . . . (21). 



Jij-j — ■ joj,-} — Jii, say . . . kL&), 



In virtue of equations (21) and (22), (19) and (20) reduce to the single equation 



(A l - A 2 ) (X + /x) /<; + 2 (X + 2/x) (C — E) = 



. (23). 



Also from (17) and (18) 



X ( r ) = u o r > e ( z ) 



w z 9 



remembering that w is odd in z and that u is not to be infinite when r = 0. 
For the stresses rr 9 zz 9 vz we find from (4), after some obvious reductions. 
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But clearly if zz is to be zero all over the plane ends we must have 



u = 






» # 



# • • » • * \ *j i /, 



k = 
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. . . (29). 



The expressions for the displacements and stresses then reduce to the following, 
writing kr = p for shortness : 
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Over the surface of the cylinder r — a, we find 
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Ay "T" -"/^ 
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Now f(z) being an even function, we can expand It in a Fourier's series between 
the limits i c ]n the form 



f( z ) —fi ) — 2a« cos 





2?i + lirz 

£iG 



1 



where 
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§ 6. Consideration of the Approximate Expressions to which the Results of the last 
Section lead, tvhen the Ratio of Diameter to Length is small. 

If we can treat the diameter of the cylinder as small compared with its length, we 
can obtain a first approximation by following the method of Professor Pochhammer 
(' Crelle/ vol. 81), and expanding A ]5 A 2 , C in powers of a, which is then a small 
quantity, provided the index n is not too large. # If we do this we find 



A- 



2 7 + 1 , b tt 2 7 - 2 

££ ....J. 1 _ 

4:j — 1 a 4fy — • 1 



A, 



a„ 



n 



27-1 , K 
47 — 1 a 4y 



9 



7 





h 



a 



u 



7 , K 7 

47 — 1 a 47 — 



r 



and, expanding I (kr) and I 1 (jfcr) in powers of r 9 and dropping all the terms except 
the first (which is really equivalent to a second approximation, since the indices go 
up two at a time), we find 



V 



u 



9( 



(4 7 - i)' 



/(c) + 2 ( a n — - 7 - b H ) cos kz 



V 



2 (47 - 1) 



f(c)+f{z)-f(c) + 



47 — 2 p 



a 






cfc 



using the Fourier expansions (38) and (39) 



r 1 
47 — 1 2 ft 



1 



(™')r=a — (2y — 1) — o (rz) rz=a dz X 2?ra 



7ra* 



Now (rz) dz X 27ra is equal to the total longitudinal pull exerted on the bar 

by all the forces on one side of the cross-section considered. It represents, in other 
words, the total tension at that cross-section. Denoting it by 7ra 3 Q, where Q is the 
mean tension at that cross-section. 



u 



r\(rr) r ^ a X , 



X + 2/x 



2{m (3\ + 2/jl) 



q 



2/i(;-$\ + 2/4 



(45), 



which shows that the radial displacement is exactly the same as if the only forces on 

a thin lamina between two cross-sections were an external radial tension (rr) r==ff 
and a uniform tension Q across the plane faces. 



* v 



For the analytical restrictions necessary in such a case, see § 28. 
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In like manner it can be shown that 



IV 



iv () z + t 



2y — 1 b n 2y \ sin kz 

""' 4y.— 1 a 4y — 1/ A; 



•'« 



^{»/w + r.(/w -/(«))*}+ j^-f ■■[*[.+(«)«* 



z 



2 7 - 1 p (rr) r 
4y - 



1 Jo 



-a 



P 



dz + 



7 



47 — 1 J fjb 



- dz 







X 



( W\= a }_+?>_ Q\ d , 

3\ + 2/A yL6 3X + 2fL /Jb J 



- f •■ 

Jo 



-(#£ 



(«), 



5^ being the stretch parallel to the axis in a cylinder which is under a tension Q 

across its plane faces and a radial tension (rr) r=a . 

Thus the longitudinal and radial displacements are, to a first approximation, the 
same as if the cylinder were supposed made up of any number of thin circular 
laminae, piled up on top of each other, the longitudinal tension in any lamina being 
uniform and giving a total tension equal to the total pull of all the external forces 
acting on the cylinder on one side of the section considered. 

Further, the shearing stress rz at a point inside is found to the same approxima- 
tion to be given by 



rz 



J) . . v 

= X (kr) ~ sin kz = %b n sin k 



z 



a fju 



so that, in the parts of the cylinder to which external shearing stress is applied, and 
in these only, there is shearing stress inside the cylinder, which shearing stress is 
proportional to the distance from the axis. 

The other stresses, rr, <£<£, zz 9 are found to the same approximation to be all 
constants for any given value of z. 



-K ° it J —f^> 



zz = 2u-S ~ cos kz 

OL 

rr = <jxj> = il/(z), 



a 



\p (z) dz = Q 



It follows from the above that the action of any radial pressure will be purely local, 
and also that, whatever the manner in which the cylinder is " gripped " and the pull 
is applied, the stress in the portions of the bar between the points of application of 
the pull reduces practically to a uniform tension. 

The above results are somewhat remarkable as tending to show how very restricted 
is the effect of local stresses, provided they leave no total resultant, and how, when 
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they do leave a total resultant, the effect of this resultant is practically independent 
of the manner in which it is applied. This is the celebrated " principle of the 
equivalence of statically equipollent loads/' which was first enunciated by de Saint- 
Venant on general physical principles, and has been considerably confirmed by 
Boussikesq's researches on. the effect of small local surface actions. 

It is to be borne in mind, of course, that the solution obtained in § 5, although 

making zz = over the flat ends, does not at the same time ensure rz = 0. In other 
words, we have a determinate system of radial shears over the flat ends, but from 
symmetry this system must be self- equilibrating. The disturbances due to it will 
therefore, by the above principle, be purely local, and, provided we remove the ends 
sufficiently far from the parts of the beam, which we desire to study, no trouble need 
arise on account of all. the conditions not being strictly satisfied. 

\ 7. Numerical Problem. Expressions for Strains and Stresses, 

Let us now return to the exact expressions and apply them to the case of a 
comparatively short cylinder. . 

Suppose that rr = all over the curved surface and that in some way, as described 

in §1, a shear rz, which we shall take uniform and equal to S, is made to act 
along two rings upon the curved surface, so that 

(rz) r=a = when — b + e < z < b — e 

z <C -~b — e 9 z > b + e 
(rz) r=:a = S when b — e < z < b + e 



We have then 



( rz) y=a = — S when — b — e < z < — b + e 



u = iv = 0, 



a n = 0, 



1 8S . 2n + lire . 2n + lirl) 

ub = -- sm. — - — -sin — . 

r (2n + 1)tt 2c 2c 

The expressions for the constants, stresses, and displacements then reduce to 



^ '-\ 



A 4S . 2n -f lire . 2n + lirb 7* 2 L -f (1 — 7) al ft \ 

At = — tt 77 sm - ■ sm L 



1 (2w + 1) TTfjL ' 2c ' 2c yct*I * - (1 + 7^) If 

A — 4S / 2 ^_+_l^ 7 27^4^1-37-6 (1 + 7) alp - (2j-_to3)Ii 

2 ~~ (2n + 1) 7T/x Sm 2^~~" Sm 2<5 7«V "" 0- + 7^) V 

C 4S . 2n + Ittc . 2^+1tt6 7*I n — 7L 

— sm — — — - sm * 



. . (47). 



h (2n + 1) w/a 2c 2c 7a 9 I 3 - (1 + T J) I, 
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rr 



V 



i 



(l+y)A 1 +(l - y )A a }I (/>) + 2 



P 



C 



A 



:p i: i^) 



cos Jcz 



r~ 



fJb 



VI 



p 



= t 



^ 



tm*4 



ZiA\.-i 



\ (p) 



1-720 

7 /<: 



Mp) 



cos kz 



&d\j 



{(2y + 1) A, - (2y - 1) A,} 1 ( p ) + r pi, (p) 



cos fc,' 



20 



( Aj + A 2 ) I x (/>) + — pl ( P ) 



?i 



h 



tlA^ip) + Tplo(/>) U ' OS ' 



. (48). 



sin h 



h 



h 



IV 



t ( A 3 I ( P ) + - pli (/>) ) &in ' 



k 



h 



i 



hi the above a is the argument of the 1-fun.ctions, unless the argument is written. 
To simplify the expressions we shall take rra = 2c, so that the length is about three 
times the radius. This makes a = 2n + 1 (n = 0, I, 2, . . . ). Further, suppose 
e = 0/6, 7; = c/2, so that the cylinder is divided into 5 zones, as shown, in fig. 2. 



'? 



+£ 



Fig. 2, 
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JS 
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The middle one from — c/3 to + c/3, unstressed ; two rings from c/3 to + 2c/3 
and — c/3 to — 2 c/3 over which a uniform shear is acting; finally, the outer rings 
2 c/3 to c and — 2 c/3 to — c, which are unstressed, Also, in order to simplify 
still more, we shall suppose Poisson's ratio to have the value 1/4, or y = 2/3. 

It may be objected, it is true, that in. many actual materials Poisson's ratio is 
not 1/4. But this is not really an objection, because the object of this investigation 
is not so much to find out the absolute values of strains and stresses in any given 
material, as to calculate the alterations in these values as deduced from the hypothesis 
of uniform, stress, and this we can best do by taking a value for Poisson's ratio 
which is, on the whole, well within the limits indicated by practical results, and which 
makes the arithmetic somewhat easier. 

If we do this and calculate the values o± the constants, we find that for the first 
10 terms they come to the following values ; 



CIRCULAR CYLINDERS UNDER CERTAIN PRACTICAL SYSTEMS OF LOAD. 163 



Table of Constants. 



/xttAj 


/XTrAr, 


//7rC 


"IS" ' 


4S ' 


m ' 


- -272602 


4- -205790 


-1 -159464 


-•142172 


— 


•0359051 


-i- -0354223 


+ -0323529 


+ 


0163035 


•— 


00534980 


+ -00492450 


+ 


•00308747 




000612345 


- -000534505 


• -- 


•000374718 
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0000532623 


- -0000285886 


— - 


0000214593 


+ 


00000237642 


- -00000413531 


— 


00000325082 


-f 


000000294829 


- • 00000162159 
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00000131798 


+ 


000000101204 


+ -000000315996 


+ 


•000000263390 


— « 


0000000175353 


+ -0000000448505 


+ -0000000381292 


— 


•00000000224045 



From these I have calculated the coefficients of the Fourier's series for the 
stresses and strains for r = 0, r = *2a, r = '4a, and r — '6a. For higher values of r 
the convergence becomes slower and the expressions more difficult to handle. In the 
case of the stresses and strains at the boundary r = a, special methods of approxima- 
tion have to be resorted to. 

The expressions for the strains and stresses are : 
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8S 



c 



fXTT" 



•00482 cos ^ + -00380 cos -^ - "00230 cos ^ 

2c 2c 2c 

•00043 cos -^ + -00006 cos -^ + . . . 

2c 2c 
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9ttz Wttz 

— -00330 cos ~ + -00084 cos -,- + -00011 cos — 



2c 



2c 



lc 



07T.i 



Ibirz 



4- -00004 cos— ;-" ;: + -00005 cos n 

2c 2c 



19 



TTi 



•0000 1 cos -— - + . . . 

ZiC 



2c 



•00002 cos 



17 



7T' 



2c 



(r = '6a) 



and in like manner for w; and the stresses. 

To save space, the coefficients of the series may be exhibited in tabular form as 
follows :— 

Y 2 
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Using the above values I have calculated the stresses and strains for the points 
r = 0, a/5, 2a/ 5, 3a/5 ? and 2 = ± (0, c/10, 2c/10, &c). These are tabulated on 
pp. 171-173. 



§ 8. Calculation of the Stresses on the Outer Surface of the Cylinder, 

Along the outer surface r = a, p = a, and we have the following expressions for 

the stresses <j><j> and zz, and the displacements u and w : rz and rr of course are 
known. 

Consider for example the stresses ( zz ) r=J « and ( <£<£ ) r==a . 
They are 



(zz) r = a= - j v-;-.T-\/-- " i 5™-; ^— L ~ ; Sin — ^~ - Sill ,. COS 



7T (ya^ * - (1 + 7a 3 )!! 2 ) (271 + 1) 



2c 



2c 



zc 
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. ■ (49), 
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4S » 27«I, 3 + 4(1 - 7) LI, - 2 (1 - 7)«I 3 . 2« + lwe . 2n + 1tt£> 2h. + Wz 

V * i _ V ' ■ U L __-..__; — — 1L oil-) -- - (sin — . P()« - 

7T 7 (7« 3 I 3 - (1 + 7« 3 ) I^) (2» + 1) 2c 2c '>- 



c 



t » 



, (50), 



Now, when a is fairly large (say > 10), I and I 2 may be replaced by their semi- 
convergent expansions : 



/ 1 / 13 13 <>2 12 <B K2 
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O , O . -j JL 



(8a) 3 3 ! (8a)» 



* « * 



/ 



(see Gray and Mathews, ' Bessel's Functions,' p. 68). From which we find that 



the coefficients of cos - — — — - in the expansions of — ( <£<£ ) r ^ a and --- (zz) r ^ a approx- 



2c 



imate to the values (remembering a = 2n + 1) : 
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Now let us write 



6 7 *V - (47 + 2) 1^ - 2 7 * V 1 



7* S I 3 - (1 +7« 8 )I 1 3 
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9, 



9 



c 



so that p f / 9 qj are comparable with the terms of the series 2 



1 



"" 35 which converge 



(2?i + iy 

fairly rapidly. We see therefore that zz and (fxf> are made up of two kinds of terms 



(a), terms of the form -— 2p« cos - — - — — and — ■ 2a,/ cos , which are abso- 

V } IT L 2C 7T * 2c 

lutely and uniformly convergent series, and (6) series, in which the coefficients are 
the approximate expressions found above. Of the series (6), those which have 
terms containing l/(2n + l) 2 or l/(2n + l) 3 are absolutely and uniformly convergent, 
This, however, is not the case with the series formed, by taking the leading terms in 
the approximation, viz. :— 
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For the series 
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(27i + 1) 2c 2c 2c 



may be broken up into the sum of four other series, thus : 
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Now it is easy to show that 
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(2/6 + 



1 / , x 
— cos (276 + 1)X 

+ 1) 
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• V 
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w 



here 



u£ 



is the numerical value of x. 
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The series on the left is divergent and log (cot — -) = 00 if x == 0. We see, 

therefore, that, at the points z = -izb ± e, i.e., wherever the shear rz changes dis~ 

continuously, the stresses zz and <jxj) become infinite. 

The meaning of this in practice would be that, as the transition from, the stressed 
to the unstressed surface becomes more abrupt, the tractions in the neighbourhood 
become dangerously large. And if the shear is applied by means of a projecting 
rim or collar of material, on which the pull is brought to bear, as in fig. 1, then 
this rim or collar must not project out of the material at a sharp angle, or in any 
way which tends to introduce a discontinuous tangential stress over the surface 
of the cylinder. This is already recognised in practice ; test pieces, which are 
thicker at the ends than in the middle, being made in. such a way that the transi- 
tion from the smaller to the larger diameter is gradual. 

The series containing l/(2n + I) 2 can also be evaluated in finite terms : 

» 1 . 2n + Ith? . 2% + lirb 2n -f Iwz 

2 --;« rv. sin ~- sin COS 

(2n + If 2c 2e 2e 






2t 
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l / Z jf~a (Jt'Z> 



16eS 
= from z ■= — c to — b ~ e 



2 
77* 



-wT (z + b + 6 ) fr° m z • = — b — e to z = — b + e 

Ibc v 



7T% 



from z = —b 4- eto z = b — e 

8e ' 



0. 



- (/; j- g _ #) fVom z ^ b — e to 2; = 6 + 

from z = b-\-e to z = c. 

Thus we have only to calculate p tt \ qj and to sum the corresponding series, the 
rest of the expressions for the stresses being reducible to finite terras. 
For y = 2/3, I find the values of p n \ qj to be given by : 
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n 




1 

2 
3 
4 

5 

6 

7 
8 
9 



Ihl 


Qn 


-•35130 


- -01184 


- -04818 


- -05960 


+ •00838 


+ -01431 


+ -00217 


+ -00436 


- -00061 


- -00135 


- -00011 


- -00025 


- -00006 


- -00015 


- -oooio 


- -00026 


4- * 00008 


+ -00024 


+ • 00006 


+ •00018 



from which the values of zz, (j><j> can he found when v = a t They are tabulated, 
with the other stresses, upon p. 171. 

§ 9. Calculation of the Displacements on the Outer Surface of the Cylinder, 
In a precisely similar manner we find for the displacements 



(«)*=* 



8Sc 



oo 



SI 



al /I x — 1 



flTT" o 



ryu^J^II-Z — 1 — ryot* 



2n + lire -'" "- 7 

— . In + kt» 2n + lirz 

2c sin. cos 



(2n + iy 



2* 



f 
^ 



2c 



8Sc'" 1 — — sin 2w + Itti sin 2wH- l7r& cos 2w 4- Iths 

/wC ^iC ZiC 



flTT 2 o 



(2rc + I) 3 



where 



u 



n 



8Sc , 2n + Itts 

'" o ^i^'U COS ~~~ 7T " 



1 ^Io/-fi — 1 \ 1 • % n + Itt^ * % n + Itt^ 

7 "" T^V/It 3 "- l"- 7*7 (2^ "+ !)" Sm 2c Sm "~2c""" 



(53) 



and is of the order 



l/(2n + 1) ; 



/ \ — 8Sc S (1 + 7) alp 3 - 7*V - 2'IJp . 2?i + 1tt0 r 2n + 1tt6 . 2n + Itt 



/wr s V (7« 3 T 3 - (1 + 7*% 8 ) (2w+ 1) J 



9 



9^ 



^ 



8Sc * 1 1 . (2n + \)ire * 2n + 1tt& . 2w 4- l7r,e 

— ; , 2 - tt ttt sm sm ~- sm 

^7r 3 7 (2^ + I) 3 2c 2c 2c 



8Sc • / l - 7 y _ 



1 . (2n + 1)« . (2n + lW? • (2?i — lVs 

sm — - — sm — sm — 



ZiC 



ZiC 



2c 



8Sc 

"i 

/AIT" 



S w« sin 



2^ + l7T^ 



» 9 



9 9 



* » 9 



» « » \ t/ rt /« 
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where 



m 



(1 + 7)*V-7«V-2I 1 T 



7 



1 



7\ 2 11 



1 



7 / a\ (2n -f I) 3 



2n -{-lire . 2n + lirh 
sm — 7 , sin - 



ZC 



ZC 



7 * 3 I 3 - (1 + 7a 2 ) I x 3 

and is of order l/(2n + I) 4 . 

It so happens that in iv the term of order lj(2n + I) 3 is evaluable in finite terms, 
and I have included it. 

It is easy to see that 



CO 



2 



1 . (2ii A- lWe . (2m + \)nb . 2n + Ittz 
: Sin sm ;- sin — 



(2» + l) f 



7T 



5 c5 



from 2: 



2c 
c to 2; 



16c 3 
2eb - £ (6 + 6 + z) 2 



7T 



1 



32c 3 



2c 
b — e 

from 2; 



2c 



& — - (3 to Z = 



fe + e 



ir^ez 



— from 2 = —h-\-e to z = b — e 



16c 3 
2<3& — i* (6 + e — 2j) 

7T 3 CZ> 



2 



7T 



32c 3 



from z = b~—e to z = b-\-e 



16c 5 



from. z = 6 + ^to^™ + c ' 



The leading series in w cannot, however, be evaluated so easily. It is seen to 
depend upon the evaluation of the series 



*x 



X 



\ log I cot - - ) dx 



. sin 2?i + lrc 

= ~ l°g I co ^ 9 ) + ^ cosec x dx. 

2 \ z / J 2 

As series of this kind are frequently turning up in investigations like the present, 

I have tabulated below the values of A — — dx and also of X ^ ^- — -- for values 

J sm a? (2% -f l) 3 

of x ranging from to tt/2 at intervals of rr/U). Intermediate values are then 
obtained by interpolation when required. 



"X 



m 



X 







1AUL 


Jii Ol ir \ -T WX — = / 

^ J sin # ^ 


\X) . 






Jj% 


/(*). 


X. 


•-238572 
•279605 
•321246 
•363596 
•406766 

! 


X. 


m. 


X. 


/(*)• 


?r/40 

2tt/40 
3tt/40 

4tt/40 
5tt/40 


•039283 j 
•078648 
•118174 : 
•157947 j 
•198050 


6tt/40 

7tt/40 

8tt/40 

9tt/40 

IOtt/40 


Htt/40 
12tt/40 
! 13tt/40 
14tt/40 
15tt/40 


•450873 
•496043 
•542417 
•590147 
•639400 


16tt/40 

17tt/40 
18tt/40 
19tt/40 
20tt/40 


•690354 
•743248 
•798291 
•855760 
•915963 
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CO 



Table of S 



sin In +ix 



o (2« + l) s 



X. 


2. 


tA/# 


2 1 


tAs B 


2. 

•78536 
•81379 
•83839 
•85938 
•87690 


eAy, 


* 


7T/40 

| 2tt/40 
3tt/40 
4-/40 
5tt/40 


•16639 

•27830 
•36959 

•44740 

•51513 

i 
I 


67r/40 
7*r/40 

8tt/40 

9tt/40 

IOti-MO 

i 


•57475 
•62754 
•67442 | 
•71602 l 

•75288 : 


IItt/40 
12tt/40 
13tt/40 
14tt/40 
15tt/40 


1 

! 16tt/40 
17tt/40 
18tt/40 
19tt/40 
20tt/40 


•89109 
•90202 
•90978 
•91442 
•91596 



We have thus the means of evaluating all those parts of the expressions which 
give rise to the most slowly convergent of the series employed. 

Taking y ™ 2/3, the values found for u n ', wj are tabulated below :— 



n. 


u n 


w n '. 





- -13933 


+ ' 03040 


1 


- -01331 


- -00634 


2 


+ -00192 


+ -00098 


3 


+ -00043 


+ -00022 


TC 


- -00011 


- • 00005 


5 


- -00002 


- -ooooi 


6 


- -00001 


- • 00000 


7 


- -00001 


- -ooooi 


8 


+ -00001 


+ -ooooi 


9 


+ -ooooi 


+ -ooooo 



Using these and the expressions given above for the finite terms, we can find the 
values of the displacements on the outer surface of the cylinder. 



§ 10. Numerical Values of the Stresses and Displacements. 

The numerical values obtained in this way are tabulated below ; I have given the 
stresses in the form of ratio (stress) /Q, where Q is the uniform tension which would 
produce a pull equal to that due to the shear S. 
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Table of Stresses. 



rr/Q. 





r = 


r — ( * 2) a r ■- 


= ( • 4) a r ■- 


= (-6) a 


r — a 


. . . . 


•22990 


•21985 j 


18587 


•11786 


•ooooo 


c/10 . . . 




22600 


•21860 ! 


19244 


•13540 


• ooooo 


2c/10 


• • 




20818 


•20842 


20227 


18098 


• ooooo 


3c/ 10 


» * f> 




17064 


•17487 ! 


18746 


21264 


• ooooo 


4c/ 10 






10688 


-10931 ! 


12179 


15054 


• ooooo 


5c/10 . 


• • 




02697 


•02239 


01531 


00700 


•ooooo 


6c/ 10 


* * 


— 


04561 


--05754 


08614 


13427 


• ooooo 


7c/10 


» • 


— 


09035 


-•10187 


13595 


18782 


• ooooo 


8c/10 


» * • 


— 


•08893 


- -09947 


12221 


14045 


• ooooo 


9c/10 


• • 


— 


•05430 


-•05936 


06827 


•06636 


•ooooo 


c 


< * • 


•00000 


•00000 


•ooooo 


00000 


•ooooo 



H/Q- 





r = 


r = ( • 2) a 


r = ( • 4) a 


r = (-6) a 


r — a 


. . . . 


•22990 


•22924 


•22568 


•21396 


•12363 


c/10 . . . 


•22600 


•22631 


•22618 


•22146 


•15364 


2c/10 . . . 


•20818 


•21209 


•21990 


•23357 


•22157 


3c/10 . . . 


•17064 


•17486 


•18827 


•21854 


•42904 


4c/10 . . . 


•10688 


•10865 


•11808 


•14158 


•25536 


5c'10 . . . 


•02697 


•02412 


•02049 


•01498 


•00162 


6c/10 . . . 


- -04561 


- -05310 


- -07023 


- -10563 


- -24868 


7c/10 . . . 


- -09035 


- -09706 


-- -11924 


- -16352 


- -41155 


8c/10 . . . 


-- -08893 


- -09653 


- -11404 


- -14443 


- -18289 


9.C/10 . . . 


- -05430 


- -05836 


- -06709 


- -07983 


- -07880 


c . . . 


•ooooo 


•ooooo 


• ooooo 


•ooooo 


•ooooo 



*/Q. 



0> 


r = 0. j r 


- (-2) a. 


r = (*4) a. 


r = (-6)c/, 


r 


= a. 


. . . . 


i 

•68906 


71895 


•81048 


•96162 


1-11724 


c/10 

2c/ 10 
3c/ 10 


* » 


• 


67272 


• 70006 


•78586 


•93696 


1 


16333 


» • * 




•63120 


65168 


•71983 


•85919 


1 


34405 


i • • 




58195 


•59404 


•63659 


•73710 


2 


02246 


4c/ 10 






53943 


•54503 


•56451 


•61686 


1 


36800 


5c/ 10 . 


• # 




50302 


•50502 


•50829 


•50813 




47865 


6c/10 . 
7c/10 
8c/10 . 






45713 


45662 


•44745 


•39955 


— 


40866 






38411 


38062 


•35965 


•27810 


__1 


05552 






27795 


•27204 


•24438 


•15199 


— 


•35747 


9c/ 10 

G 






•14532 


•14077 


•12060 


•06040 


— 


•13448 


* » 


•ooooo 


•ooooo 


•ooooo 


•ooooo 


•ooooo 
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n 



2. 


r = 0. r 


= (-2) a. 


r - (-4) cl 


r — (*6) a. 


r -= a. 


. . . . 


•ooooo 


ooooo 


•ooooo 


•ooooo 


• ooooo 


c/10 . . . 




ooooo 


02148 


•05127 


• -08883 


• ooooo 


2c/10 . . . 




ooooo 


03347 


•08205 


•15547 


•ooooo 


3c/10 . . . 




00000 ; 


03262 


•08123 


•16519 


•ooooo 


4c/10 . . . 




ooooo 


02571 


•06262 


•13159 


•95493 


5c/10 . . . 




■ ooooo ! 


02495 


•05701 


•11897 


•95493 


6c/10 . . . 




00000 ! 


03718 


•08085 


•14879 


•95493 


7c/10 . . . 




ooooo 


05812 


•12259 


•20528 


• ooooo 


8c/10 . . . 




ooooo 


07753 


•15596 


•23109 


•ooooo 


9c/10 . . . 




ooooo 


08891 


•16969 


•21918 


•ooooo 


(j • a » • 


•00000 


09231 


■17214 


•20992 


•ooooo 



S~~*. 



In the above tables it is to be remembered that rr\ zz. 



are all even functions 



c 



z 



i - , zz and (jxf> are 



of z ; rz is an odd function of z. At the points r = a, 

both = + cc 5 while at the points r = a, 2 = ± -,/, #z and ^>^> are both = ■■— oo . 

The displacements £6 and w have been compared with the corresponding total 
elongation and lateral contraction iv Q and u of the same cylinder under a uniform 
tension Q over its plane ends. 



Table of Displacements. 



W VJ . 





r = 0. 


r = (-4) a. 


r ==-• (-i)a. 


r = ( • 6) a. 


r = a 


. . . . 


•ooooo 


•ooooo 


•ooooo 


•ooooo 


•ooooo 


c/10 . . . 


•05693 


•06005 


•06988 


•08685 


•10972 


2c/10 . . . 


•11132 


•11693 


•13489 


•16752 


•22900 


3c/10 . . . 


•16235 


•16949 


•19259 


•23676 


•38253 


4c/10 . . . 


•21132 


•21915 


•24461 


•29489 


•59238 


5c/10 . . . 


•26013 


•26829 


•29467 


•34706 


•67152 


6c/10 . . . 


•30896 


•31741 


•34425 


•39538 


•68809 


7c/10 . . . 


•35493 


•36360 


•39027 


•43724 


•57421 


8c/10 . . . 


•39299 


•40162 


•42700 


•46684 


•51756 


9c/10 , . . 


•41809 


•42646 


• 45002 


•48254 


•49745 


O • . a ' . 


•42684 


•43506 


•45774 


•48725 


•49196 
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U/Uq. 



z. 






l 


r = 0. 


. 

c/10 . 
2c/10 . 
3c/ 10 

4c/ 10 

5c/10 . 
6c/10 
7 c/10 
8c/ 10 
9c/10 . 
c 






■ 


0000 

0000 

0000 

0000 

0000 . 

0000 

0000 

0000 

0000 

0000 

•oooo 



r 



1 

- ('2) a. r -- 


= (-4) a. r 


= ( * 6) a. 


r = a. 


•0449 


0375 


1341 


5789 


•0388 


0294 


1120 


5488 


•0262 


0170 


0635 


4578 


•0202 


0284 


0454 


1847 


•0314 


0856 


1206 


3709 


•0575 


1767 


2731 


4963 


•0838 


2569 


4127 ! 


5861 


•0934 


2803 


4466 


5907 


•0790 


•2313 


•3535 ! 


•3741 


•0446 


■1283 


•1880 


•1807 


•oooo 


•oooo 


•oooo 


•oooo 


i 









§ 11. Discussion of the .Results. 

The numerical results tabulated above are illustrated by the curves contained in 
Diagrams 1-6. Diagram 1 shows the radial shift, of course enormously exaggerated, 
u on the diagram, being taken as numerically equal to 2/5 ths of the radius of the 
cylinder. For convenience in plotting, the horizontal and vertical scales are not 
the same, thus a/ 5 and c/10 are represented by the same length on the diagram, 
although their actual ratio is A/tt. The same arrangement has been adhered to in 
Diagram 2. 



Diagram 1. — Distortion of a Cylinder extended by Shearing Stress applied to the Curved Surface 

(Eadial Shifts). 




(uq = 2a j '5 on the scale of the diagram.) 
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From Diagram 1 we see at once that a discontinuous change in the slope of the 
deformed outer surface of the cylinder occurs at the points z = ^ c/3, i 2c/ S, 
between which the uniform shearing stress is applied. Referring to equation (53) 
we see that at those points clu/dz changes abruptly by the value — JS//x, where 
S is the abrupt increase in the shear. This result is exhibited in the curves referred 
to, and we notice that the effect of shear, applied to the outer surface of a cylinder, 



.Diagram 2. 



-Distortion of the Cross-sections of a Cylinder under Shearing Stress applied to the 

Curved Surface. 

(w Q — c/2 on the scale of the diagram,) 




a -da -6a -4a -£& o -Ea -4a -ea -da a 

1, 2, 3, 4, 5, 6, 7, 8, 9, 10 undistorted cross-sections. 

la, 2a, 3a, ia, 5a, 6a, 7a, 8a, 9 a, 10a distorted cross-sections. 



is to depress that part of the surface towards which the shear is acting. In fact 
the greatest contraction throughout the cylinder occurs near the points z = ± 2c/3 
and appears due to this effect. 

Near the ends the cylinder broadens out again, as we should expect, though it is 
to be noted that the distorted generators meet the plane ends obliquely, which 
should not be the case if the condition of no stress over the plane ends were 
accurately fulfilled. This we know is not so : there is a system of finite shear over 

the plane ends, as is easily seen on referring to the table of rz on p. 172. This 
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system of shears is, however, self- equilibrating. The shear is zero at the centre arid 
at the circumference, and its greatest value does not exceed about 1/4 of the laterally 
applied shear. Its effects, at some distance inside the cylinder, will therefore be 
small compared with the effects of the large and unbalanced lateral distribution 
of shear. 

We notice that, for so short a bar, the lateral contraction is very much less than 
the contraction we should expect according to the " uniform tension " theory. In 
fact it never amounts to 60 per cent, of that contraction. For points deeper in the 
material, the contraction is much smaller than this. Thus, for r = (*2) a, the lateral 
contraction is 22 per cent, and for r = ('4) a it is 9 per cent, of what it should be on 
the " uniform tension" hypothesis. This seems due to the fact, in itself extremely 
remarkable, that there are considerable radial and cross-radial tensions inside the 
material. Indeed, referring to Diagram 3, we see that the radial tension amounts to 



Diagram 3 

'3Qr 



Showing Stress rr for the Cylinder under a Shearing Pull. 




about l/5th of the mean tension Q which would give the same total pull, and which has 
been consistently taken as the unit of comparison. These tensions are changed to 
pressures after passing the ring of shear, which is in accordance with the general 
compressive effect mentioned above. 

It may be noticed that the shape of the successive curves on Diagram 3 suggests 
that, as we approach the outer skin, the two bumps on either side of z = '5c would 
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■X> 



lead to infinities, or at all events, to discontinuities In the stress. In other words, 

that, though we have chosen our constants so as to make ( rr ) v=a formally zero, yet 

the limit of (rr) as given by the series Is not zero when r approaches a. This 

w r ould suggest that the series for rr, considered as a series of I-functions, behaves at 
r = a in much the same way as a discontinuous Fourier series whose general term 

Is sin nz behaves at z = tt. In fact, if we differentiate rr In the usual way with 
regard to r and then put v = a, we get a divergent series. 

It is easily seen, however, that no discontinuity really occurs except at the points 

where the shear is applied discontinuously. The general term in rr Is of the form 
(dropping irrelevant factors) : 



cos 2n + lu 
(2n + 1) 



ya[aI 1 (a)I ( / ))- / )I (a)I 1 (p)] + (l-. y )[I 1 (a)I ( /0 ) 

Up) 



a 



P 



J o ( a ) x i (p)~ 



^ 



P 



Cr a2I i( a ) — yp%(p). 



\. 



ya 2 I s (a) — (1 + yo?)I^{a) 



J 



IT 

where n = — (z i & i e )> 



Now, looking at the semi-convergent expansions for I and \ we find, putting 



/ 'Y) —I— I TTfi 

a — d =: v and p = a — 8 where S = -- — - and cl is small. 

1 2c 



I (« — S) 
X o00 



€ | 1 ~f~ It ~~ 
\ a 



1 + J- -~ 2 + terms of order S/a 3 and higher terms in S/a 



and 






e 



i + j 



2 



§ 



a 



1 



o S 



f -; + higher terms 

a" 



where e = base of Napierian system of logarithms. 

The coefficient of cos (2n + 1)^ then reduces to the form 



_1 
(2?i 4- 1) 



■f ^(a)!^)^ x — - X e s ( 1 + terms in — + terms in — -j- &c. + terms 



2 a 3 



1 



ya ^ 



i 



ill - ~f~ • . • j 

a 



+ I (a)I 1 (a)8€ 6 ( 1 + i " + terms in ~ and -^ + higher terms) > 



/ 1 x S 

+ (1 ■— y) Ii (a) I (a) {e~ 5 S ( terms of order : -- ) — — e~ s (a finite term) } 



Ii(p) 



2S 



ya a I 1 (a) < §€ s ( terms of order \ ) + -- e"' s (a finite term) 

\ a"/ a 



(ya 2 I 2 (a) - (1 + ya «)I 1 3(a)). 
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Now, if we remember that ya 2 I 2 (a) — (1 + ya 2 )I 1 2 (a) is of order yal^(a), we see 
that the successive terms in the coefficient of cos 2n+ 1 u are of the orders 

Be~ 8 Be~ 8 Be~ 8 Be" 8 Be' 8 Be~ 8 



(2n + If ' Qln + 1) ' (2n 4- If ' (2n + If ' (2n + l) 3 ' (2n + If ' 

respectively. Also in considering discontinuities, we need only consider the terms 
towards infinity, for the terms at the beginning can introduce no discontinuity. 
But clearly the series 

are of the order d multiplied by a series, which is finite and continuous up to and 
including the value d = 0. They tend therefore to the limit with d, and can 
introduce no discontinuity in the stress. 
The same will be seen to hold of the series 

Be" 8 ■"- — -- 

2 /9 ~ na cos 2n-\- 1 u, provided u ={= 0. 

If however %b = 0, we have to deal with the series 

2c (2?i + 1) 

The series under the sign of summation is divergent if d = 0. If, however, d is 
small, but still finite, the series can be summed, and we have the expression equal to 

ird. I + e- 2c 

4c' lo s [rzTt, 



e 2c 



This tends to zero when d is small, provided 

7" lo S (' """ € 2c ) j ? *^> 77 *°g 9? ^ enc i "to zero, 

which is known to be the case. Hence this series again can never introduce a dis- 
continuity in the stress. 
Now consider the series 

2 o t cos 2n + lu = — 'Sc^cos 2n + 1^. 

2% + 1 2c ' 

This is not of the same form. The series under the % is sometimes oscillatory, and 
sometimes divergent, but is never convergent, if d is put equal to zero, 
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But if d be small but still finite, the series 

ttcI and 

Jin+ivd ^— -^ (e2c — e2c')cos w 

Se 2c cos 2rc + lu ■= ~-^-~7^k px . 

1 + € H2c7 — 2e 2 ^c/cos2w 

As cZ approaches the limit 0, this series also approaches the limit 0. Hence, 
a fortiori, this series multiplied by d approaches the limit 0, and the stress is con- 
tinuous, 

(?7rdj2c 

This holds provided u dt 0. But if u = the series in question = r~~~-^ c • The 
limit of — :-— — -- when c£ = is — A. But when d = absolute zero in the series 

2 C 1 — e *tf/2c ^ 

rfiS0"~" a cos2ft+lit the series = identically. 

We have therefore in these cases a finite discontinuity in the stress. This takes 

place at the points u = 0, i.e., 2: = i ^ ± #? where the shear rz varies discontinuously. 

At all other points rr approaches the value zero continuously as we move up to the 
outer surface of the cylinder. 

Coming now to the distortion of the cross-sections, this is exhibited in Diagram 2. 
The displacements are exaggerated, as in Diagram 1, iv being taken = Jc. The 
cross-sections become hollowed out in the middle, the greatest longitudinal extension 
taking place at the sides. Another noticeable feature is that the cross-sections are 
slightly curled round the rim, except over the part of the cylinder which is subjected 

clii* du r i 7 Z 

to shear, where they slope up sharply. This follows from the fact that ••—- + t~= — . 

iVt Ltii Lb 

Thus, where rz = and da/dz > from Diagram 1, it follows that div/dr < or, 
since div/dr > nearer the centre, a maximum value of iv occurs at a comparatively 

small distance inside the " outer skin" of the cylinder. When, however, rz increases 
by S, we have seen that dufdz increases by — -| T S//x, hence dw/dr increases by f S//x, 
and is always positive at the outer surface. At the further end, where S ceases to 
act, the reverse takes place. 

It is now easy to understand why the tension is infinite at the inner end of the 
shear ring and the pressure infinite at the outer. For if we take two parallel near 
cross-sections, the one just inside the shear ring and the other just outside, the dis- 
torted cross-sections remain sensibly parallel until we approach, the outer surface, 
when they diverge sharply, if near the inside boundary, and converge sharply if at 
the outside boundary. In the one case we get an infinite extension, in the other an 

infinite compression. Hence we should expect the stresses zz and <fxf) to become 

infinite at these points, and the stress rr to vary infinitely rapidly — and this, we 
have seen, is what does actually occur. 

Further, we see that if we measure the elongation of the outer skin as is done 
with an extensometer, we shall always get too high a value for the extension. 
Referring to the table on p. 172, we have the following table of the displacements 
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Diagram 4.- 



-Showing Stress ^>c/> for the Cylinder under Shearing Pull. 
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measured by the extensometer, as compared with the displacements calculated from 
the ordinary theory, over the free length of the bar :— 



Displacements. 



Actual , . 
Calculated . . . . 
Difference , . 

Percentage correction 



(•!),. 



•10972 
•10000 
■00972 

- 8 ' 86 



(•%. 



•22900 
•20000 
•02900 
- 12-66 



,(•3) 



c. 



38253 
30000 

08253 

21-57 



(•4),. 



59238 
40000 
19238 
32-48 



(■5),. 



67152 
50000 
17152 

25-54 



We see, therefore, that in such a case very large corrections have to be applied to 
extensometer readings. 



Di 



agram 5.- 



-Showing Stress zz for the Cylinder under Shearing Pull. 
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Diagrams 3-6 give curves showing the variations of the stresses, with z, for the 
values of r equal to 0, (*2) a 9 (*4) a, (*6) a, a. I have omitted the intermediate value 
(*8) a, because the series used converge in this case inconveniently slowly, and no 
methods of approximation, such as were employed in the case r = a, are here avail- 
able. Observation of the curves for the smaller values of v will, however, in most 
cases suggest the process by which they are deformed continuously into the curve 



Diagram 6. — Showing Stress rz for the Cylinder under Shearing Pull. 
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corresponding to r = ct. In Diagram 3, of course, this is not obvious, but here, as lias 
been shown, discontinuous changes occur. In Diagram 6 it is also not quite clear 
how the curve for r= (*6)a becomes transformed into the rectangle corresponding 
to r = a. The curve for r = (*6) a has, however, already developed a double hump, 
and its righthandmost ordinate's rate of increase is fast diminishing. This suggests 
that the two humps will rise and approach each other, ultimately covering the 
rectangle, whilst the two " tails " will dwindle down to zero. 
Remarks of a similar character apply to Diagram 5. 
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§ 12. The Second Problem: Case of a Cylinder under Pressure 'whose .Ends are 

not allowed to expand. (First Method of Constraint,) 

Consider a cylinder (fig. 3) subjected to the following system of load : — 

(1.) There is no shear rz along the curved surface r = a. Over two rings of 
breadth e at the ends a radial pressure P is made to act, this pressure being so 
adjusted that there is no radial shift at the points A, B, 0, 1) ; the breadth e being 
in the limit to be made indefinitely small 



Fig. -3 



m a Q 




7Ta £ Q 



(2.) The plane ends AB, CD, are constrained to remain plane, and are subject to a 
total normal pressure rrarQ. 

The above would fit the case of a cylinder compressed between two rigid planes, 
into which shallow circular depressions had been cut, to fit the ends of the compressed, 
cylinder and prevent them from expanding. 

If we return to the expressions for the stresses in the general case, (24), (25), (26), 
and also to those for the displacements, we find that if w is to be constant for z = c 



= nrr 

\j , zzz \j^ ^2 {j 

E, ^ E 3 = E = 
( A x — A 3 ) yk + 20 = U 



6 8 9 8 



• 8 « a S * 9 \ *J *J J • 



(56). 



Also this gives rz = over the plane ends, so that we may suppose our rigid con- 
straining plane to be also smooth. 
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The condition that there is to he no shear over the curved surface now gives 



< A . + A «> + f ' W = ° 



t » • « t 



• (57), 



(writing a = ka, p == lev as before). 
From (56) and (57) 



A 1= - 



An 



C UI n 1 
Ah f 



~A 



c 



a. 



-n 



& V i, 



7 
1 

7 



y 



• • # a * t 



J 



. . (58). 



In what follows, the argument of the I-functions, when not written, will always be 
assumed to be a. 
We find 



rv = 2 (X + fx) w + \?.0 O + S 



2^0 



\ 



1 + ^) To (P) 



a 



I 



. T 

\ J i 



o + I\Ii(p) 
7 



+ P 1 l(p) 



COS / 



■*2 



Putting in this p 



- a 



(^)r=«= 2 (X + fO ^o + X^ +S 



2/aC 7a 3 y - (1 + 7a s) 1^ 



Z; 



7 aI i 



COS & 



(59), 



• * 



(60). 



Now expand the given pressure in the form 



00 



( vv) r ^ a = ; P a + 2 a w cos 



%7n 



6 



(61). 



Where a , a 1? . . . a.„, . . . are determined, P remains a free constant. 
We have at once, comparing coefficients, 



2/iC 



7 aI i 



ya?I { ? — (1 -f 7a 3 ) Ij 3 



Prt 



;/. 



* » « S * 



(62) 



2 (X -|- A 1 ) M o + ^o 



■Pa. 



* * * 



« « 



(63). 



Next we have u = when r 



a, 2 = c 







w a 



-s(-iy 







7-3 



al 



1 



T i 7 



C 

■M + '/£ a ■'■0 

tv 



Uffib 



-H 



» • 



* > 



• a 



(64), 



where 



£ = 2 ----- —— - 

1 2fxk 7« 3 1 3 • 



ii!(-Li)- 

- (1 + 7 * 3 ) I 



8 * * * 



S • 



(65). 
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This condition gives u in terms of P, and hence by (63) 



\w — P ( a Q + 2 (X + p) - 



tea 



(f'C\ 

a • • « • t U U /• 



We have now to find such a value for P that the mean pressure on the plane 
ends is Q. 

rzzdr 



ira?Q = 27r 



o 



7ra 3 (2Xu -f- (X + 2/x) w ) 



00 



+ 2ttS cos 

i « 



??,7r^ r 



{X + 2/1 }A 2 — XA- 



/> /Jo Jo 



whence, applying the well-known theorem 



(I 

dx 



(x*I H (x) ) = tf'I^x), 



we have 



■a 2 Q 



7rcr^5 = 7ra 3 (2Xw + (X 



00 

+ 2ir% cos 



%7T£ 



C 



2 ft) w ) 

(X + 2/*) A 2 — XA X - X 



20' 



" X r l + 2/.C "^ 



Z 



/• 



Using the relation 



9, 














096e$fiae 



• (67), 



and putting in for A 1; A 3 their values in terms of C, we find that the terms under 
the 2 vanish identically. Hence 



Q = ~ (2Xw + (X + 2/t) w ) 



P 



X 4- 2fA 

X 



9, 



a 



o 



f(8X + %0£ 



. (68). 



Now suppose the distribution of stress is such that rr = from % = — (c — c) to 

2; = + (c — e) and ?t = — P from z = — c to 2; = — (c — e), and from 3 = c — e 
to s=: c, we find 



^n = — e/c a, = 



//■ 



whence 



2( — 1)" . 717T6 

— sm ? 



00 

2 



c 



sin 



U7TC 



c 



alf 



~~ ' <o o ox o 

1 /a ^7r v 7a"I () - 



(1 + 7« 2 ) I, 3 



• (GO). 



When a is at all large, the terms of £ are comparable with those of the series 



c 


00 

1 


sin 


G 


M7 


n f 


->. 
"7T" 
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which is equal to 



c 



TT^fJLy 



-I 



T log (2 sin 2-) + 



'^ c x . x 7 

- cot - ax 
2 2, 



* • 



. . (70). 



(70) will give the approximate value of £ whenever ira/c is at all large. If a = 2c, 
this method of approximation will already be quite fair. 

We see therefore that if e tends to zero, £ also tends to zero, but £/e tends to 
become logarithmically infinite. 

Now from (68) 



Pa 



(X + 2/*)a + 2 ^(3X + 2yu.) 



a 



XQ 



(X + 2,0 - -^ (j) (3X + 2/*) 



Hence, since — t,/e tends to oo when e tends to zero, Pa tends to zero when e 
tends to zero. 

And similarly, for any finite value of n, Pa, u tends to zero when e tends to zero. 
But if we write down the expressions for the stresses, they are : 



ryr-/ 



Q + 



V 



7 a ~V 



Pa„7«Ii («) 

-(1 + 7^)1/ 



pi, (p) + I„ (/>) ( 2 



a. 



Ix 



7i7T# 



cos 



c 



rp/y* 



Pa + 2 



P^VyaL 



7« 3 V-(1 + 7 a 2 )I 1 2 



_(~» + i j i» W 



aT„ . 1 



I l^A + ^ + ^ 



W7T2 



COS 



y 



<£<£ — Pa + 2 - 



PtVyal! 



754 3 I 2 — (1 -f 7^ 3 ) I 1 ' 



foil 



\\ 



0^ l\Ii(f>) 



7/ P 



1 

- - lil (p) 



W7T2 



COS 



rz = 



s IXy^i 

7a 2 I 3 — (1 -f- 7a 3 ) I 1 3 



I, 



p\ fa) — y° I i fa) 



sm 



7?7T^ 



6' 



> 



. . (71). 



Now the above series are absolutely convergent for all values of r exce]3t r = a, 
where indeed they are discontinuous. Leaving the neighbourhood of r = a out of 
account, we see that for points inside the material, when the space over which the 
constraining pressure acts is indefinitely reduced, LP<x. = and 



Q 



rz 



-1/V11 



<£<£ = ; 



therefore outside the rim, where plastic deformation may be expected to occur, the 
stresses are exactly the same as on the ordinary hypothesis. 

We come then to the conclusion that this method of preventing the ends from 
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expanding is not adequate, and that to obtain any real effect, we require to make 
the constraining rim of a certain definite thickness. 

In so doing, we are really introducing an additional condition, besides the non- 
expansion of the ends, the cylinder being now, as it were, built-in. The problem as 
it stands did not appear of sufficient interest to warrant the expenditure of arith- 
metical labour upon it, so I have contented myself with stating the algebraical 
results. 



§13. The Second Problem: Constraint effected hy Shear over the Terminal 

Cross-sections. Determination of the Constants. 

Suppose now that we consider our cylinder subject to the following conditions : — - 
(i.) A total pressure rra 2 Q over the plane ends, the distribution of this pressure 

being unknown, 
(ii.) The ends constrained to remain plane, so that iv = const, when z = ± c, 
(iii.) The ends not to expand along the perimeter 

u = when r = a, z = ^ a 

This condition is satisfied by allowing a shear rz over the plane ends, its 
distribution being, however, unknown, 
(iv.) No stress across the curved surface, i.e., 

rr =: when r = a, 

rz = when r = a. 

These conditions will represent the state of things which we may expect to hold it 
the cylinder be compressed between two rigid planes which are sufficiently rough to 
prevent the expansion of the ends. 

Now, in such a case as this, it is obvious that the expressions for the stresses and 
strains as purely periodic series in z break down, for if we take the expressions (24) 

and (26) for rr and rz the condition that w = const, when z = ^ c will give us, as 
before, E = 0, and the vanishing of the stresses at the curved surface will give two 
homogeneous equations of condition between A 1? A 2 , and C. These, taken in con- 
junction with equation (35), give three linear homogeneous equations in A 1? A 3 , and 
C, which are in] general inconsistent unless A i ^ 0, A s = 0, C = 0, which would 
destroy the periodic solution altogether. 

We have therefore to assume that u and id are made up o± two parts. The 
first part, which I shall denote by U, W, consists of the periodic solution hitherto 
obtained. The second part is a finite power series in r and %. The resulting expres- 
sion is a combination of the two types of solution, which are discussed separately by 
Mr. Chree (' Camb. Phil Soc. Trans./ vol. 14). Either of these two types, taken by 
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itself, is of comparatively restricted application, but by combining the two we are 
enabled to deal with far more general problems. 
Assume therefore 

u - v r + -J, + -j- + ~ + _ + — + U (72), 

w = w^ + -^ + -* — |- — - + —r- + --— + W . . . . (73). 

O ^j ^ 4fc 

The above power series are the most general expressions of the fifth degree con- 
sistent with the conditions that u must be odd in r and even in z, and w must be odd 
in z and even in r. 

In the above we have, as before, 

U = 2j J — j}\ (kr) — - rl (&r) I cos fe * . . . . (74), 

W==s{y 8 I (*r) + f ^Ii(^)jsin^ ..... (75). 

Consider first of all the condition that w is to be constant when z— fc. This 
fixes & : 

h = 727r/c (76). 

Further we have 

F = (77), 

£D'c + P'c 3 = ......... (78). 

Now remember that u and w have to satisfy the differential equations 

, .Id du Id/ dw\ , , o \ ^ a 

( x + ^7 r &+ ^7^( r *j + ^ + 2 ^ = a 

The parts U and W we have seen already will satisfy these equations, provided 

A l — A 3 + 2C/yk = ........ (79). 

Consider therefore only the algebraic terms. Of these u r and iv {) z always 
satisfy the above equations. 
The third order terms require 

tu l (X + 2/x)+ /*D + (X + /*) D' = (80), 

2 (X + /x)D + 2/xD' + (X + 2fi)2wj = (81). 

di -B ^ 
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The fifth order terms give 

/24 \ 

(X + 2/x) [jug* + iErz 2 ) + /^Er 3 + 3/xFrs 3 + (X + p) (SEW) = 0, 

(X + /jl) (4Er% + 2F?j 3 ) + 2/xE ; # + (X + 2/i) [4i^ 3 + 3EV%] = 
which imply the four relations 

(X -p 2^£j ~ -j- fJLJh rz: Q 4 . » . . , , . . . . ( o2), 



2 (X + /*) F + 2/iE' + (X + 2ju) 4w s = ..... (83), 

4E (X + 2/a) + 3uF + 3E' (X + n) = . . . . . . (84), 

4E (X + p) + 3E' (X + 2/x) = (85). 

There is, however, a further relation to be satisfied among these constants, and 
that is obtained as follows. If we proceed to write down the expressions for rr 
and rz and to put in them r = a, we shall obtain expressions of the form 

rr = algebraic polynomial in z + series of cosines of nirz/c, 
rz = algebraic polynomial in z + series of sines of nnz/c, 

where the coefficients of cos nirz/c, sin nrrz/e, contain the two undetermined con- 
stants A { and A 2 . 

We may now proceed to expand the two polynomials in series of cosines or sines 
of nirz/c. Equating then the coefficient of each cosine and sine to zero, we can 

make rr and rz zero over the whole of the curved surface, and at the same time we 
obtain two equations for A x and A 2 . 

But it is clear that, if the Fourier expressions in the second case are to be con- 
tinuous, then the algebraic polynomial part of rz must reduce to zero when ^ = ± c, 
otherwise its expansion in sines of nirz/c is discontinuous, and at the perimeter 
of the flat ends the shear is discontinuous. This introduces infinite stresses at this 
point which render the solution inconvenient. 

Now we have at our disposal nine constants ; these have already been made to 
satisfy the seven homogeneous equations (78), (80)— (85), and therefore we are free 
to make them satisfy an eighth homogeneous equation. 

Choose then the constants so as to make (polynomial part of [du/dz + dw/dr] when 
r = a, z = ± c) zero, and we have 

Dae + Eah + Fac 3 + D f ae + EW = . . . . . (86). 

If now we express all the other constants in terms of the constant E, we find ; 
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4ryc 2 
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3 



E, 
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4 
3 



(1 +y)E : 



E' 



| r E 



and these satisfy the equations and the condition (86). 

It is noticeable that a solution can be obtained, in the form 



u 



W 



3 



Drz 2 



v + =s- + =r+u. 



ItfjS' 



3 



t(; = io Z + -~- + W, 



3 



. (87), 



which can be made to satisfy all the conditions except (86). If, however, one 
works out this solution, it is found, as we should expect, to give infinite values for 
the stresses, all round the perimeter of the plane ends. Thus, though simpler in 
form, this solution is not really simpler to work with. I have given on pp. 217-219 
the expressions for the stresses and displacements obtained from such a solution. 



§14. Determination of the Coefficients so as to Satisfy the Conditions at the 

Curved Surface. 

If we write down the expressions for the stresses, we find 



rz 



We have therefore to make 

- (D + DO az - Eah - (E ; + F) az^ 



(D + D ; ) az + Ea 3 ^ + (E' + F) at? +j+ dW 



dr 



Now we find easily 
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6' 



Hence if we expand 
S«.»sin — ■ we imd 



(D + D')az — Ea% — (F/ + F)az 3 in a series, of the form 



c 



a, 



n 



[_ (D + D') ao - Ea 3 c - (E' + F) ac 3 ] (- 1)— * x 
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The first term is zero in virtue of equation (86), and using the relations (87) 
we find 



l(mc 3 



a. 



n 



(- l Y^V^ X E X (2y+ 1), 



n 5 7r 



whence, comparing coefficients 



A, + A s + ¥ - 



(-i)""°I : (2r+i)E . . . 



. (88), 



This gives us rz consistently zero right up to the plane end, 
Ne^t we have 



( rr ) r=sa = 2 (X + fi) u + \w 
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(89). 



Now we have 
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Hence if the right-hand side of equation (89) be expanded in the form 
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we find 
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whence using equations (87) we find after long but otherwise straightforward 
reductions 
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Hence, equating coefficients on both sides of equation (89) we obtain the relations 
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[(l+y)A 1 + (l-y)AJI +2 



All C j 



(-D'Bxi^ + H*" 



W7T 



O 9 O 



« • » 



. . (91). 



and it in (91) we substitute for A 1? A 3 their values in terms of G/k deduced from (79) 
and (88), we have 



2ftC « 3 7l 3 - W -r a 8 y) 



h 



ccyl^ 



(— l) ;/ /xE X 



9. „.9 



4c% 






4c 2 



2y -f- 1 -f" ;) o (8y 4" 1) 



wwrr 



4:cl () (2y +1) 
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or 



o 



JJbQ 



(-1) 



n 






4 

(27 + l)ot + -(87 + 1) 



I x _ 4I (2 7 + 1) 



7a 3I o 2 ___ (1 _j_ ^2) Ti 2 



. (92). 



A } and A 2 being known in terms of C 5 all the constants are now determined, 
except v , iv , and E. 



§15. Determination of the Constants u , w , E. 



Let us now apply the condition that : 



f 



cW rzzdr 

J 



ITGTU, 



If we write down the expression for zz, using the expressions for Aj and A 2 in 
terms of C, we find : 



zz = 2Xt£ + (X + 2ju)w 










y W]L + (X + 2/X) — 



3E ; 



r 2 + [XD + (\ + 2 l i)w { ]z* + 2EX + — (X + 2/x) 



9 



+ — ^^r 4 + 



XF 



+ (X + 2/x) w 3 



^ 



+ 2 — - 



P-Mp) — Io(p) 



Oil, 



™ 9 



COS 



7177*2! 



C 



I (P) 



16<zc 8 7m# 



x -r-, °S r n (p) , , ioac u 

+ (2y + 1) pE 2 -^ ( - 1)" X -7-7 cos 



1 I 



W 8 7T 8 C 



2X« + (X 4- 2/a) tt' 



+ ju-E 



r 



•i 



-| C 2 + ( 2 y - 1) 



« 3 



2 



r 2 4- 



a 



I c 2 ( 2y 4 1 )] £ — 2r V J 



J- 



,i (2y-l)r 4 ' + t(2y4-lK 



00 

4- 2 



2^0 



1 /£ 



P 1 ! (p) - lo (p) 



a 



II 



2 



?17T# 



COS 



16ac 8 



+ (2.y + 1) juE S -T^.(— 1 )" X "V *~ cos 



1 Ii 



7l S 7T d c 
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r2w 

Jo 



fa. 



cl0 zzrclr =" ira % (2\u + (X + 2/x) w ) 



+ /xE X 27r 



f a 4 

ilia 



(% — 1) o t 9 "o 






'4 
.3 



(2 r + 1)6* -a 2 ] 



ft 6 2^ — 1 a'%" ~^ 



6 4 






-i 



+ aV - r ~--- : -- + 16a (2y + 1) 1 (— 1 






I ■ r J OUo 



> 



I, (a) 



J 



1 /r> Jo " " I J-i 



J J 



dp cos 



9Z7r# 



c 



But 



*« 



Ph (p) — P l o (p) 



a. 



•0 



-n 






c/i 



«%(«) -al, (a) j^O _ 2I = by (67), 

'1 J 



T /(V 1 ! J 



and 



pl (p) dp = al 3 (a), 



so that we have to make 



2X% -f (X + 2/jL)tv () 



jJLlh 



a* 



4 



<; 



2 7 






O 1 o 

a 4 + 4c- v 



a 4 2y — 1 (27 + 1) 



6 4 



o 



I 



05 / mum J yl ^fffr, 

+ 16 (2y + l)c 4 2 v ,-/ cos—" 



""\ 



(2z 2 c 3 — s 4 ) 



> = -Q. 



****€ 



Now it is easy to show that 



00 

1 



(— iy* 7i7r# 

— cdg — - 

n* c 



gV 4 



*3«.4 



IS/,4 « 0A/>2 



77T 4 



24<?» 720 



whence finally 
2Xw + (X + 2[i) w 4~ 2jaE < 



# (27 -"""* 1) 

__ 



4. i a %^ - ^ (2y 4- 1) cH = - Q . . (93). 



(90) and (93) thus give us already two equations for w , w Q , and E. We require a 
third equation. 

This is obtained from the condition that 



\U) s =a C v. 



?* as a 



VOL, CXCVTII. — A. 



2fi 
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This gives 



Da 



,9 







Ea 8 c 3 . Fac 4 






Ai 



+ * -* 1 - 







A 



9 al o 



► (-1)" = 



w a + 



"° -(1 _ y ) a? _|c 2 ] E - £ (i _ y ) E + ^ [| (1 + y) c 3 - a 2 ] E 



8 



24 



2 



+ E f 3 --Hl+y)Ea^-(2 r +l)i^E + i / i^(-l)" 



0. 



Now 



00 



t(- I) 



Z'7 



a ; 






1 



where 



£ 



oo 



, 24 



4 

(27 + 1) a + - (87 + 1) 



I] 3 -(87 + 4)1,1, 



7« 3 I 2 -L/(1 + 7 «*) 



so that £ is a known constant. 



t B 



(94), 



00 



1 



We then find, putting in for S tj its value 77^/90, 



« 







E 



1-7- r 

12 



a 



iaV+l(l+y--A-{2y + l})c* 



. (95). 



If now we write 



/= (2y -1) fi + \cW - If (2y + l)c* 

9='h(l-y-0^-i^ + j0 +y™A{2y + l})c4 . . (96), 



/t' ■«-— • ' "i "e O y 



ft- 



Q „.9,,vL 



19 ( 2 7+ 1) + -i7- 

.1. <Cl t.) 



so that y, #, A are known constants, then equations (90), (93), (95) may be re-written 

as follows : — - 

2\u + (X + 2/x) w + /x/'E = — Q ] 



t£ 



o 



+ <?E = 



2 



(X + f 1 ) w o H" ^'o + ME = 



« » « t e 



• (»?)■ 



Solving, we have : 



/xE = 



Q (27 — 1) 



(h + (2y - l)f + 2g (iy - 1) ) 



4 « « 



• • (98) : 



u 



Q.</(2 7 ~1) 







fl{h+ ( 2 y-l)f+2g(4,y-l)) 



» ♦ * « • 



. ( 9 9), 
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w = - 



Q 



h (1 — 7) +■ 2gy 



. (100). 



All the constants 
complete. 



fih + (2y-l)f+2g(4<y-l) •••••. 

are therefore absolutely determinate and the solution is 



16. Expressions for the Stresses. 



The reduced expressions for the four stresses are given below : 



rr = juE 



-r^<- 



a? 



.2 



o »0» + ^yc* f + r 2 -| (1 + y) I - (4y + 1) ~ 



+ z 2 X 2y [2c 2 - a 2 ] - r^ 1 ^) + rV (4y + 1) - 2y2 



•z- 



16ac s 



2(-l)"/IiO>) 



+ — r (2y+l)2V71 

( 7T X ' i Trl^a) 



I (p))cos 



^7T^ 



6' 



7T" i W 



( 2y + l) a + -(8 r +l)ll 1 - (8y + 4)Io 



ya 2 ll7-TT-l7" 
al n . 1 \ Ii(p) 
7/ P 



[(f + i)w»>-(f+ 



P 1 ! (/>) 



COS 



W7T2 



(101). 



^ = — Q 

+ ^-_j2 7 _-l a4+i 



3 -a 2 c 2 



if(2y+l)<r 



-|_ ^2« 



(2 7 ~ 1) 



a 2 + 



2c 3 



z2 (iK 2 y + l)c 9 ' — a s )— r 



2 7 -l 



- 2rV + §(2y -f 1) 



< 



, 16ae 3 . . » ( — 1)« I., (p) wtt2 
+ — v (2y + 1) % K - -/■ ffl cos 



7T 



\ 



. (102). 



4c% 3 7 »( — !)» ^ 

7T 2 1 ?l 2 



4 



(2y + l)a + -(8y + l) K - (8y + 4)1 



-o 



I 



T 







4* = pE - 



/ 



2 7 + 1 
12~ 



2<y 

3 



a 4 - 7 a 2 c 2 + ||yc 4 



_j_ r 



3 



(3y-l) 



— 2 



M7T.3 



COS 







4 






c 3 



vA 



57 



- (4y - 1) _- \ + 2y2 2 {2c 2 — a 2 } 






\ 



12 



+ rV(4y — 1) — 2y# 

16ac 3 



4 



7T 



(2y + 1)5, } ^—t L --rr" cos ( 

1 w° pll(#) c 



. (103). 



4c% 2 7^ ("~"1)' 
7T 3 ! n 2 



, [{(2y + l)a + ;(8y + l)]^ - (8y + 4)I 



\ 



-y« 2 (V - I1 8 ) - I1 3 

L ^+J)-w)(;-i 

2 2 



717T2 



COS 



J 
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rz 



ju,E 



{% ;(2y + l)c 3 — a % \rz + rh - |(2y + 1)? 



■t/yK 



, /.n i , \ 16ac° « ( — 1)* L(p) . ^tt5j 
+ (2y + 1) - — z~~ X -- .; ~t~- sin 



1 J 



4c% 2 7 ^ ( - 1)» 



7T" i 



4 ] 

(2y + l)a -f- - (8y + 1) > 1 L — (8y + 4)I 



w* 



7* 3 ( V - V) 



plo(p) 



I 



a 



Hp) 



rnrz 



8111 



C 



> 



(104), 



where, putting in for/, ;/, /i in (98) their values. 



E = 






27 



(4 7 _1) A ^ 



ft 



^ 



6 



a 



■1 * 



e * s 



12 



• . (105), 



§ 17. Numerical Example, 

Let us now consider a concrete example. Take a cylinder whose diameter is 
nearly equal to its length. This corresponds about to the dimensions used in practice 
for test pieces under pressure. We will take rra/c = 3 in order to simplify the 
calculation of the I-functions. 

Further, we shall assume uniconstant isotropy, so that y = '2/3. 

We then find ; 



u = 



to 







Q 



(10695),, « t 



(-090552), wj 



Q 



E= ■. (1-13842) » E 



Q - (-481906), w, 

fJLCl 



Q 



D = — ,(1-63584) 



£ = 2-036847, 



2 (1-46046), 



(1-01193), F 



m 



/xa 



t (-079057), 



3 

fJLCl 



-(1-01193), 



•-.(2-52982), 



• (106), 



Q 



D'= - ■- (1-10970), 



[Ad 



f=~ a* (-451889), g = a*("093947) , h 



a 4 (-455329) 
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and the stresses become : 






•455329a 4 - '923650r 3 a 3 + 1 '5909942V 



u 



00 



%** +\ L r s z* — |s 4 +\¥-a*2 y ~~:r- c„ cos 



1 1$ 



n 



'it 



nirz 



o 



33 + Q 



8^*( — l) tt U7TZ 

— 2-^-*« cos — 



•451889a 4 + •897748v% 2 - 2"411716^a s 



00 



(- 1)" 



r _ 2r%* + V* 4 + ¥r« 4 S >-^-/. cos 

JL ,4 1 Iv 



mn 



c 



.&&**(— 1)* ?17T3 

27 i ?r ^ £ 



<f>tj) 
fi'E 



!V1 



•455329a 4 - *359235r 2 a 2 + 1*590994^% 



/y'2/-j2 



9 



J 5 ^ 4 ~4 

-T ~ - 3* 

o 



00 ( — J } 



tt 



717T3 



A. COS 

n° o 



8^ 4 S'(— ^) n / W77 ^ 

S — 2 — 4 cos — - 

27 i % c 



r3 



/* 



E 



2'411716a%^ + r% ~~ ^rs s + \^a 4 I ( ^p, sin ^ 



00 



o ,. (— IV* . nirz 

— -A- a* 2 — -^ — q n sin — 



. (106a). 



where 



c 



Ii(p) 



w 



lo 0>) /Ii(«) 






76' 



L, \ 



^ + g- ) Il(«) " ¥ I («) 



W 



6^2 ( I() 2 _ J^) _ I3 



f — 



9» 

t 

Iln 



Io (p) I h («) 

76 



lU + ^JI^n 



P 1 ! (p) ~ h (p) ' 



al 



Hi 



6?i 3 (I 3 - 1^) - 1^ 



I 



76 
9n 



7™ + 



T _ 28T 
X I 3 X 



"■Ti W/ aI o + ^ 



Io , !\ J i (p) 






i» 



9 



_ p \ T i 



7 



1 Io (P) 



7J 



6?i 3 (I S - Ii 3 ) - I 



' 2 
1 



P» = l \{p)[h{ a ) 



2« 



7 " + 9TJ ^ - w 



pl (p) - —° I, (p) 



6n' ( V - I, 3 ) - ¥ 



/ 



(I 06&) 
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§ 18. Tables of the Constants for the special case taken. 

The values of these constants I have calculated for the values 1 to 6 of n and the 
values p = 0, a/3, 2a/3 5 a, i.e., remembering that in our case a=z Sn 9 for p = 0, 
7i, 2n, 3ft. These values are given in the following tables : — - 



Table of Constants, 



c n * 



n. 


r - 0. 


r = a/3. 


r = 2a/3. 


r = a. 


1 

2 
3 

4 

5 
6 


- -126474 

- -815103 x lO" 2 

- -485003 x 10~ 3 

- -275613 x 10- 4 

- -152381 x 10~ 5 
-•825385 x 10-* 


- -177294 

- -241966 x 10"" 1 

- -345613 x lO™ 2 

- -488500 x 10~ 3 

- -681835 x 10-~ 4 

- -943397 x lO"" 5 


- -375444 

- -144470 

- -553007 x 10" 1 

- -208129 x 10" 1 

- -776521 x 10" 2 

- -288544 x 10- 2 


- -901257 

- -929393 

- '949670 

- -961180 

- -968455 

- -973449 



&n* 



n. 


r = 0. 


r = tt/3. 


r = 2c&/3. 


r = a. 


1 

3 

4 
5 
6 


•971897 
•120292 

-117631 x 10" 1 
•945472 x 10" 3 
•678980 x 10~ 4 
•454266 x 10" 5 


1-132478 
-260646 

•586808 x 10" 1 
•115601 x 10" 1 
•208131 x 10-^ 
•353345 x lO"" 3 


1-573394 
•914018 
•534190 
•273635 
'129351 
•582489 x 10" 1 


1-960798 
1-998081 
2*319025 
2*545291 

2-702777 
2*816948 



/■»,* 
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9 



n* 



n. 


r = 0. 


r = aj 3. 


r = 2a/S. 


r - ■ a. 


1 

2 
3 
4 
5 

6 


- -787811 

- -155588 

- -176721 x 10" 1 

- -152827 x 10~2 

- -114613 x 10~3 

- -789128 x 10~ 5 


- -736099 

- -246518 

- -584824 x lO-i 

- -116088 x 10- 1 

- -209107 x 10-2 

- '354776 x 10-3 


- -324884 

- -431227 

- -326343 

- -189318 

- -963285 x 10-i 

- -455547 x 10-i 


1-638888 
2-052233 
2-400197 
2-623837 

2-774995 

2-882627 



A,,. 



»— "— — "■'" — — 

■n. 


r - 0. 


r = aj'3. 


r = 2a/3. 


v = a. 


1 

2 
3 
4 
5 
6 


•126474 

•815103 x 10-2 
•485003 x 10-3 
•275613 x 10- 4 
•152381 x 10~ 5 
•827385 x 10- 7 


•142956 

•129653 x lO- 1 
•127826 x lO- 2 
•134492 x lO- 3 
•148331 x 10- 4 
•169178 x 10~ 5 


•201175 

•397748 x 10- 1 
•991703 x lO- 2 
•275530 x lO- 2 
•814018 x lO- 3 
•250165 x lO" 3 


•333333 
•166667 
•111111 

•853333 x 10-i 
•666667 x 10-i 
•555556 x 10-i 



J n* 



n. 


r - 0. 


r = al3. 
i 


r = 2a/3. 


r = a. 


1 
2 
3 
4 
5 
6 


•971897 
•120292 

•117631 x lO" 1 
•945472 x lO- 3 
•678980 x 10- 4 
•454266 x lO- 5 


1-067168 
•179629 

•283739 x 10-i 
•414783 x lO" 2 
•586774 x lO- 3 
•817120 x 10- 4 


1-386651 
•477822 
•185745 

•707563 x lO- 1 
•266082 x 10-i 
•993188 x lO" 2 


2-042427 
1-664308 
1-682707 
1-703514 
1-716955 
1-725546 



n. 


r - 0. 


r = a/3. 


r = 2a/3. . 

V— 


r = a. 


1 
2 
3 
4 
5 
6 


• 000000 

• 000000 

• 000000 

• 000000 

• 000000 

• 000000 


•142956 

•259307 x lO-i 
•383479 x lO- 2 
•537967 x lO" 3 
•741658 x lO- 4 
•101507 x 10-4 


•402350 
•159100 

•595020 x 10-i 
•220423 x 10-i 
•814019 x lO" 2 
•300198 x 10- 2 


1 • 000000 
1 • 000000 
1 • 000000 
1 • 000000 
1 • 000000 
1 • 000000 
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9. 



ii' 







n. 


r - 0; 


1 


-•-000000 


2 


• 000000 


3 


• 000000 


4 


• 000000 


5 


• 000000 


6 


• 000000 



■r — a/3: 

■ 


r -.- 2a/3. 

i 


-•382470 . 


- -615993 


- -200638 


- -645084 


- -540924 x lO^ 1 


- -426702 


- -111882 x 10 ] 


- -230912 


- -204756 x 10-* 


- -112874 


- -349986 x 10- 8 

■ 


- -518839 x 10" 1 




• 000000 

•oooooo 
•oooooo 
•oooooo 
•oooooo 
•oooooo 



The above, when substituted in the formulae, give quite fairly rapid convergence 
when r<a, the convergency ratio being in this case less than unity by a finite 
amount. But when r = a, the series becomes comparable with the series 

2- — t — cos ;, where i is a positive integer, and. as in the first problem, a special 

procedure has to be adopted. 



§ 19. Methods of Evaluation at the Curved Boundary, 

When r = a, rr and rz are of course zero ; but the stresses zz and <£<£ require 
separate evalution. 

Now, if we use the series for I and I, in descending powers of the argument, the 
first few terms of these give a very good representation of the function when a is at 
all large, and will be quite sufficient for a>18, at which point the tables of the 
last paragraph stop. 

Replacing 1 , I T by these series, we find 



a 



nd 



1 



o 

O 



o 

O 



f n {*)~-l +^+ ^- s + 8a« + ° 



C n (a) 



l u (a) = 
e n (a) = 



9* ( a ) 



7. 

4 



7 
2' 



1 



1 



1 



o 
o 



\ 



2a 



8a 3 



3 



8a' 



ess 



165 335 

16a 9 ^ 8a-" ^ 



14 31 
+ 



a 



a" 



Ax. 

4a 3 



Jj x 



4 



49 , 157 , 503 , 



4a 



8<r 



o^a 



A. (a) = -- , #, (a) = 1 , & («) ■= 0. 

a 



\ 



* * a.v 9 



— J- 

a* 



. (107), 
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ow write 



z,/(«) 

9 J (a) 






where l/ x) (a) denotes the first four terms of the above series for 4(a), with a similar 
meaning for (j^ (a) and j^/ 4) (a). 

Then if we substitute in the equations for zz 4 Q, <f><j) we find : 



(^jLQ)-«^ 1 . 266S04a 4. 



4-41171.6 a¥ + !„¥ 



/~ 



8 



~f <X X 2 7 



- Sk ( ~~ i) w niTZ 

"i\ 2*i n GOo 



•I 



1 






G 



v^ 



r ]08^ 
i 



00 

7_ ^ 

'-4 



( - 1)" 



!0JL 

7 2 



W7T# 



» ( - If 

^ - / - cos 

1 w 4 



717T; 



ft 



COS 

c 



/> 



J 



+ — - 1 ( — if COS ~ 
1 27 i v ; 



14 



L ^ w 3 



y» fin 



n* 



4 



8a 4 

27 



00 



l) M COS 



nnrrrv. 



nir 



X 



35 



1 



ft 3 



1 1 
864 n * j. 



ow remembering that 



» (— If %7T2 



Tf 



~2 



COS 



?r 



c 



4 c 3 



7T 



12 



* - - T - cos — 



7-, 



^1 ^ 



z 2 7j- 4, 



c 4 48 n " J 1 24 



7tt 4 



720 



« (~~ If %7T# 

2 — -v~ cos — 



w 



^6 ^6 



4 

# TT 



f" lo 



yjt tj^i 



31tt 6 



c G 1440 c< 4 288 ' c 3 1440 30240 ' 



we find, using* rra/c = 3, 



(_2# _+Q)r= 

/xEa 4 



-a 



= 2-493930 



z 8 



z* 



8-201525 — + 2-251614- + '00972222 



2 l 



a" ' a' 1 ' " " a e 



O0 

4 2 ( — l) n cos 



%1XZ j j j_ 2 >' w 



c 



& 1 " 



_8_-2?! I _7 0._L. 

37 ^ +a\s 



29 16 



1 



W 5 J 
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And, in a similar manner 






•3842415 4 1-6416832 



z- 



a- 



1-1284672 



& 



a 



4 



n r rr y - 
c 
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§ 20. Calculation of the Series in the preceding Section, 
If we work out the values of//, gj, IJ we find they are as tabulated below 



n. 



1 
2 
3 
4 
5 
6 



Jn « 



+ -012367 
+ -000573 
+ '000078 
-h -000025 
+ -000009 
+ • 000005 



</«'• 




- -540511 

- -024011 

- -002538 

- -001259 

- -000217 

- -000083 


- -112667 
+ -006901 
+ -002580 
+ -000896 
+ -000381 
+ -000195 



Hence the parts of the t which depend upon /,//ft 2 , fi!/n z , gj/n % converge quite 
rapidly enough to allow us to stop after the sixth term. It therefore merely remains 
to evaluate the series 



CO / 



— l) tt TL7TZ 

—— cos — 

n° e 



and 



CO 

1 



(-!)« 



nir* 



n° 



cos 



G 



These cannot be expressed in finite terms, and although we may apply the EuleB- 
Maclatjkin sum-formula to these series directly, though in a slightly modified form, 
this sum-formula is not really of very great advantage, as its rapidity of convergence 
depends on z 9 and is such, for certain values of this variable, as to render the formula 
useless as an approximation to the remainder. As a matter of fact, however, the 
series were to be calculated only for values of z = ic/6, i being any integer from 
to 6. But for such values of z, the cosine terms repeat themselves after n = 6. 



rn 



hus. 



QO 



1 



(-1) 



tt-1 



/}■■) o 



cos 



n 



6 



1 



= cos irrlQ ( ^ ^-~~^ + ( 



m — rt 



°* L I Z^ 



1 



.m- 



m-- 



(7 + 12m) 3 , 



COS 



—p 



/ rii — ao 1 m=oo 



1 



2m)' 



:0 (8 + 12m)' 



'fll^zCO 



COS ^7T S 



1 



fii—CO 



1 



JJl: 



o (6 + 12m) : 



i + ( "~ l)! ,5o(T2 + 12m)V 



A precisely similar formula holds for 



t V f— COS "7T 

1 '/i d \ 
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Thus we see we need only work out the series 

m — oo 1 m~oo 1 

,£ (*T 12^ and „ Jo (7+ 1 2vif ' 

where s has integral values ranging from 1 to 12. 

These series are easily calculated, and, to them, the sum-formula is quite 
applicable. 

By this means it was found that 

2 t — cos( — 



n 



— cos ^ (1-000,0027 + ( — 1)'' -000,0598) 

— cos 2w/6 (-031,2519 + (-l) ; -000,0308) 
4- cos 3iv/6 (-004,1165 + ( — 1)' -000,0171) 

— cos4wr/6 (-000,9776 + ( — 1)'' -000,0102) 
+ cos 5mt/6 (-000,3207 + ( — 1)'' "000,0064) 

— costir (-000,1291 +( — 1)*' -000,0041) 



1 2>( — 1)" _1 /W\ 

and 2 ; — cos -— - 

1 n 6 \ o / 

■ = cosiV/6 (1-000,5611 + ( — 1)«' -003,1246) 

- cos2iV/6 (-125,4607 + ( — 1)«" -002,1368) 
+ cos Biir/6 (-037,4212+ (—1)'' -001,5343) 

- cos 4w/6 (-015,9496 + ( — 1)'" -001,1448) 
+ cos 5mt/6 (-008,2777 -f ( — 1)'' "000,8811) 

- cosiV (-004,8694 + (-].)' -000,6956), 

and the calculation of the stresses on the boundary then became a simple matter. 



2 d 2 
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21. Numerical Values of the Stresses, 



The values of the stresses, referred to the mean pressure as unit, are tabulated 



below : 



Table of Stresses 



-°i 



r. 



a/3 

2a/3 

5a/6 

a 



- 0. 



1 



13382 
•09971 
•00724 
•94841 
•89430 



z — c/6. 



1 • 10053 

1-00896 

-94843 

•88809 



- 2/;/6 



1-13436 I -1 



1 



10017 
01314 
95056 



- -87216 



ey In 

OC/O. 



-1 


•12145 


-1 


•08998 


- 1 


-01553 


— 


•96037 


— 


• 85845 



z -~ 45/6. 


,0 = 5c/6. 


-1-08000 

-- 1-05441 
-I -00557 

- -98465 

- -88177 


-- 1-03372 

- -96416 

- -97136 
-1-01726 
-1-04077 









= a 


| 


•68576 


' — 


-74952 


— 


•92845 


i -1 


•08211 


! -1 


-68635 







r? 



w 



r. 





a j 3 

2ft/3 



0. 



00274 
00181 
00110 
00000 



z 



c/6. 



•01414 
•01118 
•00459 
•00000 



^ 



2e/6. 



05134 
04201 
01759 
00000 



z 



3f/6, 



12416 
10220 
04481 
00000 



4c/6. 



25325 
20471 
09244 
00000 



5c/6. 



48151 
37182 
14271 
00000 



£ 



= c. 



•89668 
-65935 
•09977 
•00000 





z - 0. 


z = c/6. 


c 


z - 3c/6. 








r. 


^ = 2r;-/6. 


z = 4c/6. 


,5? — 5c:/6. 

- '48151 

- '42120 

- -28607 

- '18728 


Z = 5, 




ft/3 

2ft/' 3 

ft 


- -00274' 
•00299 
•01744 
•03265 


- -01414 

- -00737 
•01017 
•03013 


- -05134 

- -04126 

- -01461 
•02030 


- -12416 

- -10748 
-- -06480 

- =00516 


- *:io32o 

- -22288 

- -15199 

- '06344 


- -89668 

- -77728 
- -48160 

- '43801 



rz/Q. 



r. 


= 0. 


z = c/6. 


.0 = 2c/6. 


| 
z =■-= 3c/6. i £ = 4c/6. 

i 


z ~ 5<?-/6. 


i 

j 
i 

% — (\ 9 1 

i 
1 




ft/3 

2ft/3 

ft 


•ooooo 
•ooooo 
• ooooo ■ 
•ooooo 


•ooooo 

•00705 
•00425 

•ooooo 


•ooooo 

- -00388 

- -00171 

•ooooo 


•ooooo ! -ooooo 

- -02020 i - -05902 

- -02712 ! - -08645 

•ooooo 1 -ooooo 

i 


-ooooo 

- -14798 

- -19750 

•ooooo 


•ooooo 

- -35357 

- '44154 

-ooooo 
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In the above the values of zz for an additional value of r (viz., r = 5 a/6) have been 
computed in order to exhibit more clearly the variation in the pressure along the 
radius. 

The numerical results here given are shown graphically in Diagrams 7-10. We 

see at once that, save near the ends, the stresses rr, </></>, and rz do not differ very 
much from zero, which is the value they should have on the uniform pressure 
hypothesis. On the other hand, the axial pressure deviates throughout from 
uniformity over the cross-section, the total variation in any section remaining 
tolerably constant over nearly two-thirds of the length of the cylinder, and equal to 
about 25 per cent, of the mean pressure. 



Diagram 7. — Showing Stress zz for Cylinder com- 
pressed between Rough Planes (second example). 




C/6 ZC/6 3C/6 ^0/6 6C/6 

Lengths abongAxia cf cylinder 



Diagram 8. — Showing Stress rr for Cylinder com- 
pressed between Rough Planes (second example). 



wheWT= o . j 




C/6 ~~£G/6 "3C/6 ^£/& 30/6 C 

jengbhs along Axis of cy binder. 



We notice also, that, near the centre of the cylinder, the greatest pressures occur 
at the centre of the cross-section; whereas the reverse takes place at the ends, the 
pressure at the perimeter of the ends amounting to about If times the mean 
pressure. 

If we bear in mind the suggestion of the first problem of the present paper, that 
surface shear depresses those parts of the material towards which it acts, it is easy to 
see, physically, why such a distribution of pressure should be expected in practice. 
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The system of frictional shears required to prevent the ends expanding will be 
towards the centre : the parts of the material round the centre will therefore be 
depressed, and the compressing planes (supposed rigid) will have to compress 
the outer portions more than the inner, if the cross-section is to retain its original 
plane form, i.e., remain in. contact with the compressing planes throughout. It is 
thus not surprising that the greatest pressure should be at the perimeter, being in 
fact nearly 2|- times the pressure at the centre. 



Diagram 9.- — Showing Stress <j><j> for Cylinder com- 
pressed between Rough Planes (second example). 



Curve oMte wheJTr^o 



Diagram 10.— Showing Stress rz for Cylinder com- 
pressed between Rough Planes (second example). 





C/6 EC/6 3C/6 4C/6 5C/6 C 

Lengths a,Long Axis of cyiinder. 



±>05Ql 



C/6 aC/6 3C/6 4C/6 5C/6 C 

Lengths aiong Axis of cy binder. 



We see also that, near the mid-sections, the cross-radial traction changes from a 
pressure to a tension as we go towards the circumference, so that the outer parts of 
the material are subject to two pressures, parallel to the axis and the radius respec- 
tively and to a tension, at right angles to these. It should be noticed here that, in 
the diagrams, the ordinates representing the stresses increase negatively upwards. 
This has been found convenient in this case, where, owing to the general pre- 
dominance of pressures, the greater part of the stresses have the negative sign. 
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§ 22. Principal Stresses at each Point. Lines of Principal Stress. 
Now, when we have a distribution of stress 



rr, cf)(f), zz 9 rz, 
which is symmetrical about an axis, then the principal stresses are 

RR, (fxf), ZZ ; 

where <j>(j> is the same as before, and RR, ZZ are two tractions in the meridian plane, 
RR making an angle 6 with rr, where 

tan 20 = 
and the values of RR and ZZ are given by 



2zr 



. (110), 



RR 



ZZ 



rr + zz 
2 

rr -f zz 




+ ~k V ( rr *— ZZ Y + 4 (zr) J 



1 

2 




rr — zz) 2 + 4 (zr) 2 



.(111). 



Whence, using the tables in §21, we find the following values for RR and ZZ. 
compared with the mean pressure over the ends. 

Table of Principal Stresses. 









RR/Q. 








r. 


z - 0. 


z — c/6. 

- -01414 

- -01113 

- -00457 

•ooooo 


z — 2c/6. 


z = 3c/6. 


z = 4c/6. 


= 5c/6. 


# = c. 




aj 3 
2a/ S 

a 


- -00274 

- -00181 

- -00110 

•ooooo 


- -05134 

- -04199 

- -01758 

• ooooo 


- -12416 

- -10178 

- • 04405 

• ooooo 


- -25325 

- -20063 

- -08433 

•ooooo 


- -48151 

- -33690 

- -09804 
•00000 


- -89668 

- • 34800 
•09139 

•ooooo 



r. 



ZZ/Q. 





s = 0. 


z - c/6. 

-1-13436 
-1-10057 
-1-00898 
-0-88809 


z = 2c/6. 


2 = 3c/6. 

-1-12146 
-1-09039 
-1-01628 

-0-85845 


z -■=-- 4c/ 6. 

- 1 • 08000 
-1-05849 
-1-01467 

-0-88177 


z — 5c/6. 


rv O • 




-1-13382 
-1-09971 
-1-00724 
-0-89430 


-1-13322 
-1-10018 
-1-01314 

-0-87216 


-1-03372 
-0-99907 
-1-01602 
-1-04077 


-0-68576 
-1-06087 
-1-11961 
-1-68635 
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¥)' 



The values of tan 6 for the same points are given In the following table 



Table of Slope of Lines of Principal Stress (tan 6). 




z = 0. 











If 



X 



cjb 




0064 
0044 







z = 2.2/6. 


.? = 3c/6. 


z = 4c/ 6. 


&' -- OCji). 


z = c. 





















- '0038 


- "0204 


-0690 


- *2358 


- '8806 




- '001? 


- -0279 


- -0936 


- '2303 


- '4331 




















By the aid of the above, we may draw the lines of principal stress, which is done 
in Diagram 11, the slope being exaggerated in the ratio 10 : 1. In order to do so, 
we suppose the line of principal stress to remain always in the neighbourhood of the 
same generator, so that, in the above table, the values in any row apply to the same 

Diagram 11. —Lines of principal Stress for Cylinder compressed between Eough Planes (second example) 















J8 o 



Q 






OS 



(drawn bo base osc). 
■Una ofsbresa for r= a&fa 
fgmmi bo b& aao'scf). 




C/6 2>G/6 3C/6 -4-C/6 SCfs C 

Lengths abonO Amis of cylinder. 



& 






line of principal stress. This of course is not correct near the ends, but it is sufficient 
for our purpose, the diagram being merely intended to show the general course of the 
lines of stress. These are sensibly parallel to the generators throughout the middle 
part of the cylinder, but slope outwards near the ends. 

Diagram 12 # shows the distribution of the lines of equal principal stress ZZ 

* In Diagram 12, as in several others, a is represented as having the same value as c, so that the horizontal 
and vertical scales are not strictly the same, but differ in the ratio ir/3. This has been done for con- 
venience in plotting. 
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inside the cylinder. I have chosen the stress r LZ in preference to the others, as, 
except in the neighbourhood of the centres of the plane ends, it is everywhere the 

Diagram 12. — Distribution of Principal Stress ZZ inside the Cylinder (case of compression between 

Rough Rigid Planes). 




-$0-65^07570 
lines are drawn at intervals of 05 Q of the stress. The critical line corresponds to ZZ = 1*01 Q nearly 



greatest of the three principal stresses, and therefore ZZ will be the important 
quantity when we come to discuss the maximum stress. 

The diagram has been constructed by careful interpolation from the values 
tabulated above, and from it we see that the surfaces of equal principal stress in the 
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cylinder are in general made up of three sheets, and they fall into two classes : 

(a) those for which ZZ has a value less than a certain critical value, which, as nearly 

as I can find out from graphical methods, is about 1*01 Q, and (b) those for which ZZ 
has a value greater than l'Ol Q. 

The surfaces (a) consist of two solid caps or buttons, round the centres of the end 
sections, together with a hollow cylindrical shell surrounding the middle of the 
cylinder. For values sensibly < '9 Q the latter sheet disappears, and only the caps 

remain, their volume gradually dwindling down to zero as ZZ falls to *686 Q. 

The surfaces (b) consist of an elongated core, resembling a cylinder closed by 
curved ends, surrounding the centre of the compressed block, together with two 
annuli at the ends, as shown in the figure referred to. 

The critical surface ZZ 1*01 Q consists of two nearly plane sheets, roughly co- 
inciding with the cross-sections Z = ± 5 c/6, and one cylindrical sheet, which bends 
inwards towards the end, though without completely closing in, and which roughly 
coincides with the cylinder r = 2a/3 over the greater part of its surface. 

§■23. Application to Rupture. Distribution of Maximum Stress, Strain, and 

Stress Difference, 

In considering what happens when a material breaks, we have to ask, first of all, 
whether it be brittle or ductile. In the first case, the law of stress to strain will be 
approximately linear up to the point where rupture takes place ; in the second case, 
the stress-strain relation remains approximately linear until a point is reached (called 
the yield-point) at which a large and sudden change occurs in the stress-strain curve, 
after which the material becomes sensibly plastic, so that rupture finally takes place 
after a large permanent deformation. 

In applying an elastic theory to practice, we can, in strictness, treat of rupture 
only in the case of a brittle solid. Even then it has to be borne in mind that the 
mathematical theory of strains — upon which the equations of elasticity depend— 
requires the strains to be so small that their squares are negligible. It is possible 
that, even in the case of the most brittle solids known, this condition may cease to 
hold before rupture occurs, although the stress-strain relation may continue to be 
linear, Nevertheless, the calculated values of even the breaking strains in a material 
like cast iron, for instance, are so small as to render this unlikely. 

For a ductile metal, such as mild steel, the elastic results only tell you where the 
material will begin to take permanent set. 

In the case of stone or cement, however, to which the present results would be 
applied, there seems to be no definite yield-point or elastic limit, the material being, 
in fact, only imperfectly elastic throughout. Still, we may consider that the results 
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of elastic theory give in such a case an indication of the state of stress when the 
specimen breaks. 

There are three distinct theories, both of rupture and of failure of elasticity. 
According to Lame and Navier, failure occurs when the greatest stress at any point 
exceeds a certain limiting value. This is also often taken as the criterion of absolute 
rupture. According to Saint- Venant, the maximum strain, and not the maximum 
stress, is that which determines failure. Finally, a theory has lately been put forward 
by various elasticians to the effect that failure occurs when the greatest shear at any 
point, that is, the greatest principal stress-difference, exceeds a definite amount. 

Professor Perry has proposed another criterion, suggested by the angle at which 
compressed cylinders shear (see 'Applied Mechanics/ pp. 344-345), namely, that 
rupture takes place when s — \xfp exceeds a certain value, where s is the shear across 
any element of area at a point, p is the normal pressure on this element of area, and 
jx is a constant depending on the material. This theory, however, need not concern 
us so much, as it appears more specially applicable to the final breakdown of ductile 
materials long after they have become plastic. On the other hand, it has been 
shown by Mr. J. J. Guest (' Phil. Mag./ July, 1900) that the beginning of plasticity 
was very probably determined by the maximum stress-difference. 

Let us now proceed to apply these three criteria, namely, those of the maximum 
stress, maximum strain, and maximum stress-difference to the cylinder in the present 
example, and see what results they give us, on the hypothesis that for materials like 
stone and cement, plastic yielding and rupture are simultaneous. 

Consider first the greatest stress theory. This would make failure of elasticity 
first begin to occur round the perimeter of the plane ends, and that as soon as 
1*68635 Q > a certain limiting value S . If the pressure be uniformly applied, and 
the ends expand, we get failure of elasticity when 

Q > S . 

Hence the apparent strength of a cylinder tested in this way would be about '593 
of the strength of a cylinder tested under a distribution such as is usually assumed. 

Further, if we consider the regions where the stress is greater than a given value S, 
we find that they consist of separate spaces, which join on to each other as S 
diminishes, the critical value for which this occurs being given by S = 1*01 Q, The 
regions of greatest stress consist therefore of a central core, which spreads out into a 
sort of hollow cone near the ends. If then we suppose fracture to occur over regions 
of greatest stress, we see why it is that the material breaks off in conical pieces at 
the ends. 

Consider now the greatest strain theory. Let T } , T 2 , T 3 be the three principal 
stresses, and s v s 2 , 5 3 ^ e corresponding stretches. 

Then , i== ^ Ti __>_ (Ti+T> + t 8 ) 

2 E 2 
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so that the greatest s will correspond to the greatest T, if T b T 3 , T 3 have the same 

sign. This is our case everywhere, except in cases where <f*<f> > 0, and then <fxf> is so 

small that it still leaves the strain corresponding to ZZ numerically the greatest. 
We have then, remembering we have assumed X = /x, to investigate the values of 

ZZ - 1 ( ZZ + RR + ft ) = 2/xs, . 
This will be proportional to the greatest strain, except near % = ± c, r = 0, where 

TIE - |( ZZ + RR + #) = 2 /^- 

should be taken. It is found, however, tliat at this point the strain is comparatively 
small, and the maximum strain there is a matter of indifference. 



Table of s g /s, where s = maximum stretch under the same uniform pressure. 



r. 


# = 0, 





1-13245 


a/3 


1 • 10000 


2a/3 


1-01133 


a 


-90246 



Z = c/6. 


£ = 2c/6. 


z - 3c/6. 


z = 4c/6, 


1-12729 

1-09589 

1-01035 

•89563 


1-10755 
1-07935 
1-00509 

•87724 


1-05938 
1-03756 

•98813 
•85716 


•95338 
•94751 
•94446 
•86591 



5c/6. 



77296 
76590 
86416 
99395 



z 



c. 



•23743 

•39034 

•78311 

1*57685 



It we take therefore the " greatest stretch " theory, failure of elasticity still occurs 

at the perimeter of the ends, but this time only when the stress is ------- (limiting 

r ,/ 1-07/ x & 

stress in the case of uniformly compressed cylinder), so that although the apparent 
strength is less than in the uniform case, it is greater than if we adopt the " greatest 
stress " theory. 

The lines of equal principal stretch s ff /s are shown in Diagram 13. They are 
drawn for only one quarter of the meridian plane, the rest being symmetrical. They 

present the same general characteristics as the curves of equal stress ZZ, with this 
difference, that the critical line corresponds to s, = 4 915 s\ Again, the caps or buttons 
at the ends are far larger ; so that if pieces are cut out, they will be considerably 
larger than on the "greatest stress" theory. Also, looking at the inclination of the 
lines joining the corners to the critical points, i.e., the points where the two branches 
of the critical line intersect, we see that the fragments, if approximately conical near 
their base, will probably be cut off at a much higher angle than in the previous case, 
Let us now proceed to consider what happens if we adopt the third or u greatest 

stress-difference " theory of rupture. It is easy to see from the tables of RR, ZZ, 

and </><£ that the greatest stress-difference is either RR-ZZ or <f><f> -ZZ. In the 
sixteen cases tabulated, for which z > c/8, the first of these is the greatest stress- 
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difference, and in the twelve remaining cases the second is the greatest, although, 
as a matter of fact, the two stress-differences, for these twelve cases, do not diverge 
very much. 



Diagram 13. 



-Distribution of Principal Stretch, s z , inside the Cylinder (case of Compression between 

Rough Rigid Planes). 




^r 



The number corresponding to each line =• the value of s z /$ for that line. 

- - , critical line. s z /s -■= *9 15. 

The actual greatest stress-difference is given in the following table :-- 

Table of (Maximum Stress-difference) /Q. 



r. 


z = 0. 


• z ==" c/6. 


z = 2c/6. 

1-08188 

1-05892 

•99853 

•89246 


z = 3c/6. 


z = 4c/6. 


z = 5c/6. 


5! = C. 




a/3 
2a/3 

a 


1-13108 

1-10270 

1-02486 

•92695 


1-12022 

1-09320 

1-01915 

•91821 


•99730 
•98861 

•97223 

•85845 


•82675 
•85786 
•93034 

•88177 


•55221 

•66217 

•91798 

1-04077 


- -21092 

•71287 

1-21101 

1-68635 



Here again plastic deformation will first occur round the perimeter of the ends 
when Q = (;7T Q ^r~ ) of the value it should have, on the same theory of rupture, in 

order to produce failure of elasticity in a uniformly compressed cylinder. 

So far, then, this theory leads to the same results as the maximum stress theory. 

Diagram 14 shows the distribution of maximum stress-difference. The lines of 
equal maximum stress -difference present very similar characteristics to those of equal 
maximum stretch. The critical line corresponds to a maximum stress-difference = 
•933 Q. 

Remarks similar to those used in the last case apply in this, 
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Diagram 14. — Distribution of Principal Stress-difference inside the Cylinder (case of Compression between 

Rough Rigid Planes). 




The number corresponding to each line denotes the value of the ratio (maximum stress-difference : Q) 

for that line. 
„.__..--„_ - ? critical line. Stress-difference =-= *933 Q. 



§ 24. Distorted Shape of the Curved Surface. 
If we work out the values of u when r = a, we find, after some reductions. 



U r ~a 
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Now in the particular case we are dealing with, if a be large, r„ (a) approximates to 
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Write then 
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r »( a ) = 2-i;+'i£' + S + ? ''» • • • • ■ • ( 114 )> 
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and substitute in (112). We find, putting in for 2 -cos---- ■* and 2 -- cos — : their 
known values. 
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where, putting in now ira = 3c, we have the following values for ?\/ : — 



n. 



£i 



3 



4 



5 



(J 



■r,,,' (a). . 

-■42430 
-•01807 
- -00223 
-•00063 
- -00024 
-•00010 



Whence, using the methods of § 20, we find the following values for u/u , where 
u = lateral expansion of a cylinder of the same dimensions under the same total 
pressure uniformly distributed : — 





'U/U(,. 


V/ • » » e 


. -97861 


c/6 . 


. -96325 


2c/6 . . 


. . -91118 


3c/6 . . 


. . -80402 


4c/G . 


. -61209 


00/ O • . . 


. -30863 


G • * 


. '00000 



These results are exhibited in Diagram 15. We see that the cylinder has a single 
bulge in the centre, and indeed it is easy to verify that {dhi/dz 2 \ <0 by differentiating 
equation (115) and putting in the numerical values. 

Diagram 15. — Showing Distortion of Curved Surface (second example). 
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This agrees with the figures shown by Bach in his ' Elasticitat und Festigkeit ' 
(figs. 2 and 3, Par. 11) of the distorted shapes of such cylinders. There are, 
however, in the possession of Professor Karl Pearson, at University College, 
London, specimens of iron which have been strongly compressed, so that the strain 
has been large and permanent and the meridian section of the distorted curved 
surface is an obviously wavy curve, with two maxima, one on either side of the mid- 
section. With regard to the apparent disagreement here between theory and practice, 
I would observe that these specimens have been subjected to enormous stresses, for 
which the equations of elasticity certainly do not apply, and probably are not even 
an approximation ; in the second place, the specimens I have seen are longer, 
compared with their diameter, than the cylinder of the present numerical example, 
so that it is easy to see why the conclusions above need not apply to these specimens. 

§ 25. Apparent Young's Modulus and Poissons Ratio. 
We find that the total shortening of the bar 
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Hence the apparent Young's modulus E Y ' = Qc/(w)- r 
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Now if a be small, i.e., if the bar be very long 
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retaining only terms of 4th order, 
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where i> is the greatest integral value for which a = vna/c makes the above approxi- 



mation suihcient 
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Hence if a/c is very small, v may be large, and thus although the first term in 
a 4 £/6 cancels -4-5 c 4 , yet the terms in curled brackets will become indefinitely large 
compared with the other terms. Thus, for a very long cylinder, 






/x,(47 - 1) __ ™ 



7 



Y? 



where E Y is the true Young's modulus. 

On the other hand, for a very short cylinder, c/a is very small, and £ being of the 

order — r -7; X 2 J s , the two leading terms in the numerator and denominator 
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of the right-hand side of (117) are negligible and 
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This is identical with the modulus of compression for a cylinder which is prevented 
from expanding laterally by a constant pressure applied to the sides. So that we see 
that for a very flat disc the effect on the modulus of compression is the same, whether 
the lateral expansion be prevented by means of shearing-stress over the flat ends, or 
by hydrostatic pressure over the curved surface. 

The apparent Young's modulus for intermediate cases will be between those two 

values (these, for uniconstant isotropy, being E Y and 6/5 E r ). Thus, in the given 

example, where y — 2/3, 

E/= l-0498E y . 

Poisson's ratio comes out to be apparently '2690 instead of '2500. 

Thus the errors in the values of Young's modulus and Poisson's ratio, as deduced 
from an experiment with cylinders under the given conditions, will be 5 per cent, 
and 7 *6 per cent, respectively. 



§ 26. Solution involving Discontinuities at the Perimeter of the Plane Ends. 

In § 13 it was stated that a solution, obtained by methods strictly analogous to 
those used in that section, but which neglected the condition that the shear rz should 
be continuous at the perimeter of the plane ends, could be found. 

It seems of interest to give, for purposes of comparison, the expressions for the 
displacements and stresses, deduced from this solution. They are : 
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The method by which the constants are obtained is precisely the same as that 
*ed in §§ 14, 15. 
The expressions for the stresses are given below : 
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It is now easy to see, if we bear in mind that when n and therefore a is large, CI 
remains finite, as appears on examination of (121), how it is that the stresses at the 
boundary become infinite. 

For in both zz + Q and <j>(f) we have, when r = a, terms of order 1/n, when n is 
large. These are of alternate sign, C containing ( — l) n . But if z = ± c they 
become all of the same sign, and the series become logarithmically infinite, 

§ 27. Summary of Results, 

Looking back upon the results obtained, we notice : 

(a.) That the three solutions we have been considering successively are only the 
simplest of an infinite series of solutions, which are continually growing more com- 
plicated ; for we need not necessarily stop, as has been done, at terms of the fifth 
degree, but might go on to terms of any degree in r and z and thus construct, as it 
were, solutions of successive orders. We should then have an infinite number of 
free constants, which might be determined by introducing further limitations at the 
plane ends, such as, for instance, restricting n to be zero at every point and not 
merely along the perimeter. 

The analytical complexity of such a complete solution would, however, be very 
great, and would render it quite beyond the reach of arithmetical expression, and 
consequently valueless for the purposes of the engineer and the physicist. No 
attempt has therefore been made to develop this solution, although, as an analytical 
possibility, it appears interesting. 

(6.) That the different solutions all agree in giving the perimeter of the plane 
ends as the locus of the points where the elastic limit will first be passed, one of 
these solutions actually making the stress infinite at this perimeter. 

In the more important solution, however, where continuity and finiteness are pre- 
served, the conclusion still holds, and, further, is independent of whatever theory of 
tendency to rupture we adopt, whether we suppose it due to maximum stress, to 
maximum stretch or squeeze, or to maximum shear or stress-difference. 

(c.) That in the numerical example considered, plastic deformation begins to occur 
round the perimeter for a stress between 2/3 and 1/2 of that which is required to 
cause a cylinder under uniform pressure to pass the elastic limit. 

This is apparently in contradiction with the results of engineering experience, both 
Unwin and Perry stating that blocks of stone or cement, pressed between millboard, 
which hinders the expansion of the ends, show greater strength than the same blocks 
when the ends are allowed to expand. 

The key to this appears to be found in a remark of Unwin, which Professor Ewing 
confirms, that the lead sheets do not merely allow the expansion of the block, they 
force it, i.e., lead in its plastic state will exjDand more than the stone or cement would 
do laterally under a uniform axial pressure. 

2 F 2 
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But the solution, when the ends are compelled to expand by a given quantity a , 
is easily deducible from that given for non-expanding ends. Thus, let u u u\ be the 
values of u and iv in the case worked out, when Q = 1, and let us write 



Xr 



U = P (^-7^— TTTT) + B^u 






Then this will satisfy all the conditions, provided 

P _\_ _ 

/, 6\ + 4/1 X ^ " a ° 

and P + R = Q. 

Therefore « r . /,-. H^i^x . , \\ 

i/ = — + (Q - ^° £ (6X + 4/x) ) 7/ l9 

giving the solution under a given mean pressure Q, which produces a flow a of a 
lead plate, and thereby constrains the ends of the test piece to expand by that 
amount. 

The principal stress at the perimeter of the section is now 

P + (1-686) R = 1*686 Q - -686 P, 



and if P, i.e., a , be made large enough, this can be made much smaller than Q. It 
begins to be smaller than Q as soon as the expansion induced by the flow of the lead 
is greater than the natural expansion of the stone under uniform pressure. 
On the other hand, the principal stress at the centre of the plane ends is 

P + '686 R = '686 Q + -314 P, 

and this again, may be made great by making a large. 
The principal stress- differences are : 

at the perimeter 1/686 Q — '686 P 

at the centre P - -211 R = - -211 Q + 1-211 P. 

Hence we see that whatever theory of failure we adopt, if the ends are forced to 
expand, so that P > Q, the material first becomes plastic (or else breaks) at the centre 
of the cross-section, the strength of the test-piece diminishing as P increases, but 
having no definite value. That some such thing as this does really occur in practice 
is very well shown by the results published by Unwin (' The Testing of Materials of 
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Construct ion'), where blocks show less strength when three sheets of lead are 
introduced between the compressing planes and the test piece than when one sheet 
only is introduced, the lateral flow being greater in the first case owing to the larger 
amount of lead. 

It would seem, therefore, as if the true strength of a cylinder were really greater 
than its strength as tested either between millboards or between lead sheets, and 
not, as Professor Perry states in his ' Applied Mechanics/ equal to the strength 
shown in the lead test — this test, as we see, leading to results that are not definite, 
but vary with the expansion of the lead. The millboard test, however, which is 
advocated by Unwtn, should give a constant value, although it is not the value which 
would hold for a cylinder under uniform pressure. 

(d.) Diagrams 12-14 suggest an explanation of the fact that, when short cylinders 
are strongly compressed between very hard surfaces, pieces are sometimes cut out at 
the ends of an approximately conical shape. The same occurs when spherical pieces 
of metal, such as ball-bearings, are compressed between parallel plates. This is usually 
explained by saying that the material breaks along the planes of principal shear. On 
the other hand, it may be argued simply that rupture should take place over the 
regions of greatest stress. These are near the perimeter at the ends, and gradually 
close in upon the centre, forming hollow caps. 

Further, in the case of the lead tests, where P > Q, this state of things is 
reversed, and the material should give way from the inside, so that we should expect 
it to split axially, and possibly along meridian planes as well. That this is what 
really occurs can be verified by referring to the figures in. the chapter on testing of 
stone in Un win's ' Testing of Materials of Construction/ 

(e.) The results both of this and of the first problem show us how unreliable any 
experiments on short cylinders must be, which have in view the determination, by 
tensile strain, of either Young's modulus or Poisson's ratio. Thus any results 
obtained in such a case without the dimensions and the mode of application of the 
stress being exactly specified, would not justify us in general in drawing any 
conclusions as to whether a given material possesses or not uniconstant isotropy. 

§ 28. The Third Problem. Case of Torsion. Expressions for the 

Displacement and Stresses, 

I now proceed to consider a case where u and tv are zero, that is, where we have 
to deal with the solution in v, which we have seen is independent of the others. 
We have in the notation of § 3 

(£3 + 1)2) v — o . 

Hence, excluding K-functions, since the solution must be finite and continuous at 

the origin, we have 

v = £ (A;,, sin kz + B, ; cos kz) I } (kr) , 
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Now, if the torsion be symmetrical on either side of the mid-section, we have 
v = 0, when z : = 0, therefore 

B w = . 

But also 7 s dv u . 

02; ~jLtv = 0at the ends. 

Therefore cos he = or he = 2n+lw/2 . 

Iso 

r0 = ur -- - = 2uA ;/ sm A;; . r - — — — - 

' ' dr \ r ' ' dr v 



= S/xA w sin /^ . H 2 (&r) 



By the well-known property of the I-functions 
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Now suppose that the cylinder is subjected, to a certain system of transverse 
surface shears, so that rcfy can be expanded in a Fotjriek's series in the form 
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T^ ^ *1 (P) 2^1 •+ 1.7T2' , 

</>£ = S c„ r— - cos — ^;---. ....... (1-29), 







h (*) 2, 



where p a have the same meanings as before, viz., fa% &a. 

In the case where a/c is very small, or the cylinder is very long relatively to its 
diameter, we may obtain a first approximation by retaining only the first terms in the 
expression for the I's, Proceeding as in § 6, we find 
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we see that 
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r<£ = r 2 t/; (#) /cr, 
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Now if M be the total torsion moment up to any cross-section. 



M = 27ra 2 i/*(z)c£ 
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= torsion at the point = r. 
— 2M 



Therefore 



rra 



M = fi X — X r and <£# = /xrr . 



But these are the actual formula) connecting the torsion with the applied couple 
and with the shear across the section for a circular cylinder. 

We see, then, that the usual formula) continue to hold, to the first approximation, 
when the forces applied to the surface of the cylinder vary with z 9 provided we define 
our torsion-couple at any section (much as the bending moment at any section of a 
beam is defined), as the couple of all the external applied forces to one side of that 
cross-section. 

It is interesting to note also that, to this approximation, there is no distortion of 
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the cross-sections, v/r being constant for the section, Straight radii therefore remain 
straight radii. 

Further, r<j> = (r 2 /cfi) (its value at the boundary). 

In other words, the transverse shearing-stress across cylinders coaxial with the 
given one is zero for sections where there is no such external applied stress, and 
for other sections diminishes rapidly along the radius as we go inwards, so that near 

the centre it is always small compared with <f)Z. 

There is one very important point to be noted with regard to this method of 
approximation : p and a increase with n y and therefore, however small a/c may be, so 
long as it remains finite, we still reach a value of n> for which it is not justifiable to 
take for I x and I 3 the first terms of their expansions In positive integral powers of the 
argument. If, however, we stop at the vt\\ term, where v is finite, then if R, is the 
remainder after v terms of the series 
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and if, on the other hand, the numerical value of the difference 7™^ — v^j-z < € for all 
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values of n not greater than v, where e is a quantity which depends upon c/a, and 
which can be made as small as we please by making e/a small enough, then the 
fence 
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denotes the numerical value or modulus of x, and .11/ is the remainder after 



v terms of the series 



Z c n 4r 2n + 1 nrz 

2 zr -v- COS - 
k d 4 



,v 



imlO 



Now if both the original series and the approximate series are uniformly 
convergent, then by giving v a certain value, | R„ | and | li/ 1 can both be made less 
than a certain small quantity 17/3 which tends to zero when v tends to infinity, and 
that for all values of z. 

Now make c/a so small that e < t]/3p, which we can always do so long as v is finite. 
Then the difference between the two series is numerically < 17, and the approximation 
holds. 

If, however, for any value of z y it becomes impossible to assign an upper limit to 
11^ or it/, i.e., if either series cease to be uniformly convergent, then we should have 
to increase v indefinitely in order to make jll„l < 77/3, and therefore to modify e, so 
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that no limiting value of c/a (which should, of- course, be independent of z) could be 
found, and the approximation need not necessarily hold. As a matter of fact, it is 
shown in § 29 to fail for particular cases. This is true d fortiori , if either series cease 
to be convergent at all. 

The same remarks apply in their entirety to the process of approximation given in 
§ 6, and further, to the approximate expressions given by Professor Pochhammer in 
his investigation on the bending of beams ('Crelle/ vol. 81). 

§ 29, Special Case of Two Discontinuous Rings of Shear. 

Suppose that we have the following system of values for <f)Z : — 

(f)Z = T if c — e O <e, 

(f>z=- it — a -f e <z < c — e, 

(f)Z = — T if — c <z < — c + e, 

so that we have a cylinder twisted by two equal and opposite rings of transverse 

shear extending over lengths e of the cylinder, near the ends. Then we find 

easily 

_ 4T . . (2n + l)ire 

(2n -f l)7r x 2c 

with the following values of the displacements and stresses : 

V ~f vVn + lf^~~ L > I,(*) Sm 2c Sm 2c 
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<f)Z = Z 7^— 7-t— (— 1) — - sm -~ cos 



. . ( Idl |. 



o (2^+ 1)tt v 7 I 2 (a) 2c ^ 2c 

Now it is easy to see that in this case the conditions for uniform convergency are 

satisfied, except at the boundary, and except with regard to the stress r6, whose 
approximate expression is not uniformly convergent, being in fact discontinuous for 
2 = ± (c — e). 

At the boundary, I x (a)/I 2 (a) tends to unity with n, it's approximate expression, 
when a is large, being 

I-i(«) -. . 3 - 15 , 15 / x 

_iA.../ — 1 J J_ 4_ 4„ /l qo\ 

I 8 («) ~~ ^ 2« ^ 8a 3 ^ 8« 8 -r ....... . U^/- 

Hence i? is always uniformly convergent and its approximate expression likewise, so 
for it the approximation, for sufficiently small values of c/a, holds throughout. 
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For r<f> the series is non -uniformly convergent for r ~ a in the neighbourhood of the 

sections z = ± (c — e) owing to the series having a finite discontinuity. For ^ the 
approximation certainly fails, for r = a, in the neighbourhood of z =:- ^ (c — e) for 

part of the expression for ^ is the series % l/(2n + 1), which is divergent. 



Hence, it we are in such a case to use our approximations for the stresses, we must 
exclude the sections where the applied stress is discontinuous and their immediate 
neighbourhood from consideration. 

It will be found that similar remarks apply to the example of pull given in § 7, 
and also to the example given by Professor Pochhammbr in his paper on bending (loc. 
cit.) y in which he also deals with discontinuous systems of stress, so that his approxi- 
mate expressions leave us in the dark as to what does really happen at points of 
support, the cross-sections in the neighbourhood of such points being, for reasons 
analogous to those developed above, excluded from the regions where his approxima- 
tions hold. 

Before proceeding to an actual numerical concrete case, we may notice that <j>z 
becomes infinite at the points z = ± (c — e) for the causes stated above. Hence any 
discontinuity in a system of transverse shears applied to the surface of a cylinder, or 
any such shear transmitted by a grip applied to a portion of the material projecting 
at a sharp angle, will produce in the neighbourhood a very great stress across the 
section. It would seem, therefore, that a cylinder treated in this way would be likely 
to experience plastic flow, or to rupture, not in the middle, but near the points where 
it is seized. 



§ 30. Approximations on the Boundary when the Cylinder is short. 

When the cylinder is short, so that a becomes rapidly large, we may use the 
method employed in §§ 8, 9, and 19, availing ourselves of the approximation (132). 
We then find :— ~ 



^_4TC J_7T/ _ g- 



e 



log. 



cot 



7T 

¥ 



z — e\ 



e 



«Jm 



4Te 



r f^ 1 "'"?) 



Vl-i+- e )log, 



4 



c 



cot 



7T 



1 
2 



f 



x cosec x 



dx — |r x cosec x dx 

~ J n 



1 



Z + 6 



+ f 



rn 
17T 



fJLG 



ez (from z 
and 



to % = g — e) 



ez 



|- {z — c + e T (from % = a — e to z — c) 



CIRCULAR CYLINDERS UNDER CERTAIN PRACTICAL SYSTEMS OF LOAD. 227 



T^ » / 1 

+ 15-^S(-1)"' 



flTT o 



MMIM 



(2w + l) 4 ' (2?i4-l) 5 y 



271 + 17TC . 2% -f- llTZ 

sm 1 sm - 



2c 



2c 



" 1 " ^tt 3 f (27i + 1)2 [ I 8 («) 2a 8«* 8« s 



15 15 1 . 2n + lire . 2u-\-1ttz ^ ^^\ 

sm — ___ — - sm # # (133), 



2c 



9/* 



4> 



T 
-log, 

7T 



tan 



/ 



7T 7T 3 — e 



4 4 c 



tan ( I 



7TZ + e 

1 ~~ 



+ s 



ttT 



4 



C 



<3 (from 2 ■ = to & = c -— e) 

and 

c — 2; (from ^ = c — e to z = c) 



T 00 

7T o 



n 



1 1 \ . (2n +• V)m 2% + Itt? 

X^TTTy 5 + <>T+ i)v sm ""IT " cos 2c 



3 



4T- (-!)» /I^)^ _ 
"*" ir o (2n + 1) \I»(*) 2a 



15 

8 a 2 8a 8 



15\ . 27i -f Ittc 2n 4- Itt^ ^ ._ A , 

sm cos — — (134). 



2c 



2c 



where in the last t in both equations only a comparatively small number of terms 
need 5 be taken. 



§31. Numerical Example. Values of the Coefficients and of the Displacement 

and Stresses. 



Take a cylinder such that wa/2c — 1, and let e = c/2, so that the distribution ot 
stress is as shown in fig. 4. Then a = 2n + 1, and we find : 



v 



4v/2Tc 



fJL7T 



2 



(v sin itz/2c — v l sin 3nz/2c — ^ 2 sin 5irz/2c + v 3 sin 7irzj2c + t; 4 sin 9irz/2c 



^ 2 /2T 

r ^ = 1__ (f. Q s i n ^^c — t x sin 3irz/2c — . . .) 



Vgsin 11ttz/2c — . . .) 



<^£ ■=: ~ ( o * COS 7T2:/2C — ^ COS 37T^/2c — . . .) 

the law of the signs being obvious, and the v's, t% and s's being given below : 
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From these, using the formulae of approximation given in the last section, when 
r = a, we obtain the values of v and the two stresses, I have tabulated them in the 
form v/v and (stress) /T where v 0> T are the greatest values of the displacement and 
of the shear respectively in a cylinder of the same length, subject to a uniform 
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torsion over its whole length, the total couple applied being the same as in the 
present case. We find T = yura = ?rT and v = rca = ttcT/jjl for the given example. 

Tables are given on page 229. 

Looking at these tables we see that, over the length free from external applied 
shear, the strains and stresses inside the cylinder are sensibly the same as what they 
would be on the hypothesis of a uniform torsion. Outside zjc = *5 the torsion couple 
diminishes, and the stresses diminish in consequence. 

It is interesting to compare these results with those that we should have obtained 
if we had supposed the approximate results given on p. 223 to hold good in this 

case. Denoting by v\ r<j>\ <f>z' the values of the displacement and stresses calculated 
on this hypothesis, we have — 
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5 32. Discussion of the Results. 

From the above tables we see that the radii in each cross-section do not remain 
straight lines, but assume distorted shapes, which are shown, on a very exaggerated 
scale, on Diagram 16, where, for each of the ten cross-sections, the curve ol 
(v — v')/v Q9 which indicates the deviation from the straight line in the distorted 
form of the radii, has been plotted. The variation from the straight line increases 
rapidly as we approach the region where the stress is applied, as can be seen from 
curves (l)-(4) on Diagram 16. On the other hand, towards the ends, the distortion 
remains fairly constant. The distorted radii meet the bounding circles at right 

angles when r<f> = at the surface, but they meet it at a finite angle where 

r<f> = T. 

From the values of <j>z and r<j> we see that as soon as we get at all away from the 

ends the conditions that r<j> = 0, <j>z = \xjtz, v = rrz, which hold for uniform torsion, 
are very closely satisfied, and that, more generally, except where the abrupt change 
takes place in the shearing stress at the surface, the approximate expressions given 

in § 28 do not differ widely from the true expressions, the law that r<j> varies as 
the square of the radius being, near the ends, tolerably well verified. It is to be 

noted also that, where the approximations would give a discontinuity in r<f> inside the 
material (viz. at z = '5c), the true values are almost exactly the mean of the two 
discontinuous values obtained from the approximate formulae assumed correct. 

In like manner cj>z is nearly the same as <f>z\ except near z = '5c, where, as we have 
seen, an infinite stress really occurs, of which the approximations give no hint. 

We note, however, that <j>z does not strictly vary as r all over the section, being 
smaller than should be expected inside and larger at the boundary. 

The theoretical result, that the stress is infinite where the transverse applied shear 
is discontinuous, throws much light on the case of a cylinder whose cross-section 
abruptly changes, as in fig. 1, with the difference that now the stress applied to the 
collar is transverse. We see that in such a case we should expect the material to 
give way at the points of sudden change. This conclusion is in accordance with 
practical experience, the tail ends of propeller shafts, for instance, breaking almost 
invariably in this manner. 

§ 33. General Conclusion. 

This example concludes the series of three which it was proposed to treat of The 
object has been to obtain a clear idea of the effects of certain surface distributions of 
stress which come much nearer to the cases arising in practice than does the uniform 
distribution ordinarily taken. 
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Diagram 16.— -Showing Distortion of a Radius originally Straight in the case of Torsion produced by 

applying Shearing Stress to the Curved Surface* 
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The curve corresponding to the section zjc = n/10 is numbered n. The first four curves have had the 
ordinates exaggerated in the ratio of 10 : 1. They are shown by the dotted lines, and to them refer the 
numbers in brackets. 
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No doubt the cases treated involve somewhat arbitrary conditions, not strictly 
obtained in practice, but it appeared useful to ascertain how far they gave results 
diverging from those which would be found on the ordinary hypothesis of uniform 
extension or torsion. 

This furnishes us with a test of how far we may accept de Saint- Venant's 
principle of " equipollent " systems of load for a bar whose length is gradually made 
smaller compared with its diameter. The results we have here obtained indicate 
that, as we go away from the points of application of the stress, a " uniform " solution 
is reached much sooner in the case of torsion than in that of either tension or 
pressure. 

With regard to the arithmetic of the paper, the results have been as far as possible 
checked. It is believed that they are correct to the number of figures given, but 
owing to the slow convergence of certain of the series, accumulated errors may in 
some cases affect the last and even the second last figure. Even this, however, would 
not sensibly disturb the conclusions. 

For the I-functions the tables in Gray and Mathew's " Bessel's Functions " were 
used, but the range of the tables is so limited that a large number of these functions 
had to be independently calculated. The semi- convergent expansions were employed, 
the argument being large in each case. 

My very best thanks are due to Professor Ewing for his unfailing kindness in 
coming to my aid with suggestions and advice. 
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